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ABSTRACT 

Prompt gamma-ray burst (GRB) emission requires some mechanism to dissipate an ul- 
trarelativistic jet. Internal shocks or some form of electromagnetic dissipation are candidate 
mechanisms. Any mechanism needs to answer basic questions, such as what is the origin 
of variability, what radius does dissipation occur at, and how does efficient prompt emission 
occur. These mechanisms also need to be consistent with how ultrarelativistic jets form and 
stay baryon pure despite turbulence and electromagnetic reconnection near the compact object 
and despite stellar entrainment within the collapsar model. We use the latest magnetohydro- 
dynamical models of ultrarelativistic jets to explore some of these questions in the context 
of electromagnetic dissipation due to the slow collisional and fast collisionless reconnection 
mechanisms, as often associated with Sweet-Parker and Petschek reconnection, respectively. 
For a highly magnetized ultrarelativistic jet and typical collapsar parameters, we find that 
significant electromagnetic dissipation may be avoided until it proceeds catastrophically near 
the jet photosphere at large radii (r ~ 10 13 -10 14 cm), by which the jet obtains a high Lorentz 
factor (y ~ 100-1000), has a luminosity of Lj ~ 10 50 — 10 51 ergs -1 , has observer variability 
timescales of order Is (ranging from 0.001-10s), achieves yOj ~ 10-20 (for opening half- 
angle 9j) and so is able to produce jet breaks, and has comparable energy available for both 
prompt and afterglow emission. A range of model parameters are investigated and simplified 
scaling laws are derived. This reconnection switch mechanism allows for highly efficient con- 
version of electromagnetic energy into prompt emission and associates the observed prompt 
GRB pulse temporal structure with dissipation timescales of some number of reconnecting 
current sheets embedded in the jet. We hope this work helps motivate the development of 
self-consistent radiative compressible relativistic reconnection models. 

Key words: accretion discs, black hole physics, galaxies: jets, gamma rays: bursts, MHD, 
instabilities, relativity, methods: numerical 



1 INTRODUCTION 

Gamma-ray bursts (GRBs) are thought to originate from core- 
collapse events or compact object mergers leading to magnetars or 
accreting black holes capable of launching ultrarelativistic jets. The 
prompt emission from standard cosmological long-duration GRBs 
has an energy of about 10 51 erg over a few seconds that is beamed 
into a jet with an opening half-angle of a few degrees (Frail et al. 
2001; Bloom et al. 2003). The prompt emission is typically pre- 
sumed to occur in internal shocks (Sari & Piran 1997). The inter- 
nal shock model is reasonable because such shocks are expected 
in an unsteady outflow, which then has an observed variability 
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on timescales related to the central engine. Within the collapsar 
(Woosley 1993; Paczynski 1998; MacFadyen & Woosley 1999) or 
other GRB models, the observed variability may arise indirectly 
from activity near the central compact object, indirectly from en- 
trainment driving propagation instabilities (Aloy et al. 2000; Zhang 
et al. 2004; Morsony et al. 2007; Wang et al. 2008; Bucciantini et al. 
2008a), or directly from relativistic turbulence (Lyutikov & Bland- 
ford 2003; Narayan & Kumar 2009; Zhang et al. 2009), 

However, the internal shock model has some unresolved prob- 
lems. For example, highly relativistic relative motion between in- 
teracting shells is required in order to efficiently generate photons 
(Kobayashi & Sari 2001; Maxham & Zhang 2009) ; only a small 
fraction of electrons should be accelerated in order to obtain con- 
sistency with the observed peak energy (Shen & Zhang 2009) ; in- 
ternal shocks produce a steeper spectral slope than observed (Ghis- 
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ellini et al. 2000; Asano & Terasawa 2009) ; the afterglow energy 
would dominate the prompt energy as opposite to observed (Will- 
ingale et al. 2007) ; particle-in-cell (PIC) simulations show that 
shocks are dominated in energy by Maxwellian electrons whose 
emission may be inconsistent with the observed double power-law 
Band function (Giannios & Spitkovsky 2009) ; and if GRB jets 
contain strong toroidal fields, then the shocks should inefficiently 
dissipate the kinetic energy (Kennel & Coroniti 1984; Zhang & 
Kobayashi 2005; Narayan et al. 2011) and inefficiently accelerate 
particles to high energies (Sironi & Spitkovsky 2009). One alterna- 
tive to the internal shock model includes dissipation from nuclear 
and Coulomb collisions (Beloborodov 2010a). 

In light of these issues with shock models and since GRB 
jets are expected to be highly magnetized, it is interesting to ex- 
plore alternatives for which electromagnetic dissipation directly 
leads to acceleration and emission (Romanova & Lovelace 1992; 
Thompson 1994; Usov 1994; Levinson & van Putten 1997; Lyu- 
tikov & Blackman 2001; Lyutikov 2006b; Thompson 2006). Such 
investigations are motivated by studies of pulsar winds (Kennel 
& Coroniti 1984; Coroniti 1990; Lyubarsky & Kirk 2001; Kirk 
& Skjasraasen 2003a; Nagata et al. 2008), field reversals in active 
galactic nuclei (AGN) jets (Lovelace et al. 1994, 1997), and mag- 
netars (Thompson & Duncan 1995; Lyutikov 2003, 2006a). Elec- 
tromagnetic dissipation mechanisms are appealing because they do 
not require the generation of highly relativistic relative motion be- 
tween shells of matter or highly relativistic turbulent motion. Dissi- 
pation can instead proceed in situ. Indeed, high-resolution 3D mag- 
netohydrodynamical (MHD) simulations of jet propagation show 
that toroidal magnetic fields shield jets against many of the shear in- 
stabilities noticed in purely hydrodynamical simulations (Keppens 
et al. 2009; Mignone et al. 2010), which suggests that it is diffi- 
cult for a magnetized jet with a strong toroidal field to produce nu- 
merous shells moving at varying relativistic speeds as required by 
the internal shock model to generate temporal variability and effi- 
cient dissipation. Further, modern simulations of magnetized GRB 
jets suggest that even if relatively relativistic shells were generated, 
internal shocks would be unable to explain the high efficiency of 
prompt GRB emission (Narayan et al. 2011). This implies that lo- 
cal electromagnetic dissipation may be the only efficient dissipation 
mechanism possible in highly-magnetized jets. 

Further, several observations point to a requirement of strong 
electromagnetic fields in GRB jets. The typical absence of a ther- 
mal photospheric emission component requires the jet to contain a 
significant low-entropy (e.g. magnetized) component, and the pres- 
ence of non-pair-producing GeV photons requires emission to be 
at larger radii than expected in the internal shock model (Zhang & 
Pe'er2009). The typical non-detection by Fermi of a GeV spectrum 
excess due to synchrotron self-Comptonization (SSC) in prompt 
GRBs can be achieved by a highly magnetized jet since fewer elec- 
trons are required to support a strong electromagnetic field, as com- 
pared to the number of electrons required by the internal shock 
model that should generate an SSC component (Fan 2009). Pos- 
sible measurements of highly polarized gamma-rays from GRBs 
could require an ordered magnetic field in the emitting region (Lyu- 
tikov et al. 2003). The absence of bright optical flashes in the very 
early afterglow may require the reverse shock region to be some- 
what magnetized (Mizuno et al. 2009; Mimica et al. 2009). 

While electromagnetic fields are commonly understood to 
be dynamically important in the jets produced by GRB engines 
(Narayan et al. 1992; Thompson 1994; Usov 1994; Vlahakis & 
Konigl 2003; Lyutikov 2006b; Uzdensky & MacFadyen 2006, 
2007), only recently have self-consistent magnetohydrodynamic 



(MHD) simulations been able to model the formation and large- 
scale axisymmetric structure of such jets (Komissarov et al. 2009; 
Lyubarsky 2010a), including what happens beyond a stellar en- 
velope in the collapsar model (Tchekhovskoy et al. 2010b). Such 
large-scale idealized MHD jet simulations are complemented by 
small-scale general relativistic MHD (GRMHD) axisymmetric and 
3D simulations that consider the role of magneto-rotational insta- 
bility (MRI) driven turbulence (Balbus & Hawley 1991) and mag- 
netic field geometries close to the black hole (McKinney 2006b; 
Beckwith et al. 2008; McKinney & Blandford 2009) and by re- 
alistic GRMHD simulations of the engine in the collapsar model 
(Mizuno et al. 2004; Barkov & Komissarov 2008; Nagataki 2009). 
Of relevance to GRBs is that a dipolar field appears required to 
launch a jet from the black hole (Beckwith et al. 2008; McKinney 
& Blandford 2009), yet magnetic field advection from a presuper- 
nova core collapsing onto a black hole leads to too small of a mag- 
netic flux to generate the power of a cosmological GRB (Komis- 
sarov & Barkov 2009). The MRI must then be invoked (Akiyama 
et al. 2003), but that only generates small-scale (not dipolar) field. 
One solution is for a magnetic dynamo to generate quasi-periodic 
large-scale dipolar fields over tens of dynamical times as occurs for 
toroidal fields in accretion disks (Davis et al. 2010), but this neces- 
sarily implies one must understand electromagnetic dissipation at 
the interface between each successive flip of the jet's dipolar field. 

However, it remains difficult to understand how the electro- 
magnetic field dissipates in GRB jets. The most advanced simula- 
tions of relativistic jets from collapsing stars forming black hole ac- 
cretion flows use the MHD approximation (Komissarov & Barkov 
2009), so they cannot incorporate all the plasma effects observed 
in PIC simulations that are required to self-consistently study mag- 
netic reconnection. Even the most advanced PIC simulations are in- 
applicable to GRB jets that require relativistic reconnection in the 
presence of pair creation and annihilation, neutrinos, photons, ra- 
diative cooling in both optically thick and thin regimes, relativistic 
compression, etc. There is still no coherent physical theory of rel- 
ativistic reconnection (see, e.g., Blackman & Field 1994; Lyutikov 
& Uzdensky 2003; Lyubarsky 2005; Jaroschek & Hoshino 2009; 
and also section 5.5 of Uzdensky & MacFadyen 2006). Therefore, 
quasi-analytical approaches must be used to investigate the role of 
electromagnetic dissipation in GRB jets. 

In prior quasi-analytical work, electromagnetic dissipation 
through magnetic reconnection in GRB jets has been suggested as 
a possible source of acceleration and emission at large radii from 
the central compact object (Thompson 1994; Spruit et al. 2001; 
Drenkhahn & Spruit 2002a; Lyutikov 2006b; Giannios & Spruit 
2006). It was considered problematic for ideal MHD processes to 
efficiently accelerate and collimate jets, so non-ideal electromag- 
netic dissipation was favored to accelerate GRB jets (Thompson 
1994; Drenkhahn & Spruit 2002a) or be a source of energy for 
emission in GRB jets (Thompson 1994). Prior works also assume 
that the flow is already relativistic (y » 1) and continues to speed- 
up with radius as based upon one-dimensional jet models. The dis- 
sipation rate was assumed to be Alfvenic or some fraction of it, 
which assumes a fast reconnection mode. 

However, it is now known that a relativistic quasi-conical out- 
flow (i.e. split-monopole near compact object and otherwise un- 
confined) accelerates efficiently near the polar axis in ideal MHD 
(Tchekhovskoy et al. 2009), and any sufficiently globally colli- 
mated outflow is efficiently accelerated (Komissarov et al. 2009; 
Tchekhovskoy et al. 2010b). This dramatically changes the behav- 
ior of the Lorentz factor as a function of radius compared to prior 
one-dimensional models (for a review, see Spruit 2010.) In addi- 
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Figure 1. Reconnection switch concept: Collapsar model or some other sys- 
tem produces a jet (with opening half-angle 8j) corresponding to a gener- 
alized stripped wind containing many field reversals that develop into dis- 
sipative current sheets. The jet collimates due to confinement by the stellar 
envelope out to a radius r mono , after which the jet (having become ultrarel- 
ativistic) cannot expand laterally and becomes a nearly radial (monopolar) 
flow. The figure shows toroidal field polarity reversals produced by dynamo 
processes near the central engine or by turbulent entrainment at boundary 
layers. In the reconnection switch model, these current sheets avoid signifi- 
cant dissipation while in the collisional regime until they reach sufficiently 
large radii where the plasma becomes collisionless and fast collisionless re- 
connection is triggered. The curved dotted lines denote that the striped wind 
continues down to the central engine. This paper focuses on the dissipation 
processes within current sheets that ultimately lead to prompt emission. 



tion, jets following ideal MHD acceleration near the compact object 
start with y « 1 until reaching the Alfven surface. So, if dissipation 
were always nearly Alfvenic, then (unabated by time dilation) fast 
reconnection could occur near the jet base. Reconnection near the 
jet base can severely quench the Blandford-Znajek effect (Beckwith 
et al. 2008; McKinney & Blandford 2009), which may be required 
to drive the relativistic electromagnetic jet in the first place. Also, 
if reconnection were quite close to Alfvenic, then the jet could dis- 
sipate deep inside the photosphere leading to a hot fireball instead 
of an electromagnetic jet. This would lead to a dominant thermal 
spectrum, which is inconsistent with observations. 

Prior quasi-analytical works have also investigated the poten- 
tial role of the current-driven magnetic kink instability, which is as- 
sumed to operate rather efficiently on Alfvenic timescales (Spruit 
et al. 2001; Giannios & Spruit 2006). However, modern stabil- 
ity analyses have shown that the kink instability operates some- 



what slower than comoving Alfven timescales for collimating flows 
(Narayan et al. 2009), and much of a relativistic quasi-conical flow 
can be out of causal contact with itself to eliminate the instabil- 
ity except within narrow regions right around the polar axis that 
contain little electromagnetic energy flux (Spruit et al. 2001; McK- 
inney & Blandford 2009). 

Prior work also considered other aspects of electromagnetic 
dissipation. For example, a violation of ideal MHD occurs when 
the plasma density drops below the Goldreich- Julian density, and 
this leads to dissipation (Lyutikov & Blackman 2001). However, for 
reasonable lower-limits on the baryon-loading (Levinson & Eichler 
2003), such a violation of ideal MHD only occurs at r » 10 19 cm 
(Spruit et al. 2001). Reconnecting layers in accelerated jets have 
been found to be potentially unstable to magnetic interchange in- 
stabilities leading to enhanced reconnection (Lyubarsky 2010b), al- 
though they considered a sharp pressure boundary instead of a dis- 
tributed radiative photosphere as present in GRB jets. Turbulence- 
driven electromagnetic dissipation from magnetic instabilities in- 
duced by internal shocks has also been recently considered (Zhang 
& Yan 2010), although the conditions for turbulence in highly mag- 
netized jets remains undetermined. 

One aspect of electromagnetic dissipation not investigated 
so far is the role of critical transitions between collisional and 
collisionless reconnection and its association with, respectively, 
rather slow resistive-MHD reconnection (perhaps as slow as Sweet- 
Parker reconnection; Sweet 1958; Parker 1957) and fast Petschek 
or Petschek-like reconnection (Petschek 1964). (For a review, see 
Priest & Forbes 2000.) 

Prior simulations and laboratory experiments show that colli- 
sional plasmas avoid the fast Petschek type reconnection in favor 
of the slow Sweet-Parker type reconnection. For collisional plas- 
mas the Petschek reconnection configuration promptly collapses to 
a thin Sweet-Parker layer (Uzdensky & Kulsrud 1998, 2000). No 
fast (inflow velocity of order the Alfven velocity) Petschek recon- 
nection is seen to occur in collisional plasmas (Biskamp 1986a; Ma 
& Bhattacharjee 1996; Ji et al. 1998; Uzdensky & Kulsrud 2000; 
Breslau & Jardin 2003; Malyshkin et al. 2005; Kulsrud 2001). 

On the other hand, a collisionless plasma allows for non-ideal 
MHD effects that can force a plasma into a fast Petschek-like re- 
connection regime occurring on the ion skin depth scale and op- 
erating independently from the classical resistivity. (For a review, 
see Kirk et al. 1994; Kulsrud 2005; Zweibel & Yamada 2009.) 
Simulations of pure pair plasmas also show fast reconnection but 
on the electron skin depth (instead of the ion skin depth) as per- 
haps due to the electron pressure effect (Bessho & Bhattachar- 
jee 2007). Petschek reconnection plays an important role in the 
Sun, where nanoflares remain the most plausible source of coro- 
nal heating (Parker 1988a; Klimchuk et al. 2009). In fact, the solar 
corona may even exist in a balanced marginally-collisionless state 
governed by the density-controlled transitions between these two 
modes of reconnection (Uzdensky 2007c, a). There is also observa- 
tional evidence that the transitions between the collisional and col- 
lisionless reconnection regimes may also control coronal heating 
in black-hole accretion disks (Goodman & Uzdensky 2008). Solar 
and stellar flares could also be triggered by these critical transitions 
(Cassak et al. 2006, 2008). Therefore, it is important to determine 
whether a reconnecting system is collisional or collisionless in or- 
der to determine whether current sheets operate in the fast recon- 
nection mode (Malyshkin et al. 2005; Cassak et al. 2005; Yamada 
et al. 2006; Uzdensky 2007c,a; Bessho & Bhattacharjee 2007). 

There do exist purely resistive-MHD alternatives to fast 
collisionless Petschek reconnection, including due to externally 
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driven MHD turbulence (Lazarian & Vishniac 1999; Kowal et al. 
2009b; Loureiro et al. 2009) and plasmoid-dominated reconnec- 
tion (Loureiro et al. 2007; Kowal et al. 2009a; Samtaney et al. 
2009; Cassak et al. 2009; Huang & Bhattacharjee 2010; Uzden- 
sky et al. 2010). However, relativistic Poynting-flux dominated jets 
are strongly electromagnetically-dominated and relativistically ex- 
panding flows, and hence they may only be susceptible to turbu- 
lence or similar cascades in regions that remain in causal contact 
across the jet - such as very close to the polar axis where there is 
little total electromagnetic energy and so little possible emission. 
These alternatives require further consideration in such relativistic 
regimes. 

The present paper builds upon three novel key elements. The 
first key element of this work is that we consider the role of transi- 
tions between slow collisional and fast collisionless reconnection. 
The fact that the reconnection rate depends upon the collisional- 
ity of the plasma is particularly important for GRB jets that are 
collisional at the jet base and out to large radii. We compute the 
properties of the collisional and collisionless reconnecting layers to 
determine which reconnection mode dominates at each radius and 
angle within the GRB jet. Once collisionless reconnection is ini- 
tiated at some "transition radius," reconnection proceeds at a rate 
much faster than the collisional rate and can initiate prompt GRB 
dissipation and emission. 

The second key element of this work is that the collisional 
layer's properties are determined as due to collisions among species 
of ions, photons, electrons, positrons, and neutrinos. The current 
layer structure is treated similarly as done for radiative accretion 
disks (see, e.g., Kohri et al. 2005). We also make use of relativis- 
tic reconnection work by Lyubarsky (2005), which shows that the 
relativistic Sweet-Parker layer has a reconnection rate as expected 
from non-relativistic theory. We incorporate radiative effects into 
the reconnection switch model by following Uzdensky & McKin- 
ney (2011), who showed that the current layer compresses in the 
strong cooling limit of non-relativistic radiative reconnection. 

The third key element of this work is that the large-scale jet 
is modelled using the latest ideal MHD models of ultrarelativistic 
jets rather than assuming inefficiently accelerated one-dimensional 
flows. We also extend the study of Spruit et al. (2001) and deter- 
mine a range of possible ways for field polarity reversals to occur 
that lead to current sheets in the substructure of the jet. 

A basic version of the overall argument of the paper is pro- 
vided in §2. The GRB jet structure is described in §3, some possi- 
ble field substructures are described in §4, the reconnection models 
are described in §5, results for GRBs and other jet systems are de- 
scribed in §6, a discussion is in §7, and conclusions are provided in 
§8. In Appendix A, the generalized full jet structure solution that 
is used throughout the paper is obtained. In Appendix B, the equa- 
tion of state for all species within the radiative current layer is pre- 
sented. In Appendix C, the collisional and collisionless reconnec- 
tion models are discussed. We assume a flat space-time in spherical 
polar coordinates (r, 6, <p), an orthonormal basis for all vectors, and 
Gaussian-cgs-Kelvin-radian units. 



2 BASIC ARGUMENT 

In this section, a basic argument is presented to demonstrate the 
existence of a reconnection switch mechanism that leads to dissi- 
pation near the GRB jet photosphere. 

Figure 1 shows a basic picture of the reconnection switch 
model for a jet containing multiple current sheets each correspond- 



ing to a layer wherein oppositely-directed magnetic field lines are 
dissipated. Near the central engine, dissipation proceeds via slow 
collisional reconnection due to the high collisional rate. This al- 
lows the electromagnetic field to avoid significant dissipation de- 
spite the presence of current sheets, and this allows the generation 
of a baryon-pure strong electromagnetic field for launching an ul- 
trarelativistic jet. At large radii, pairs annihilate and the densities 
go down leading to infrequent collisions. Collisionless plasma ef- 
fects (e.g. due to electron-proton or electron-positron decoupling) 
can then initiate the much faster collisionless reconnection mode 
(operating in a Petschek-like geometry) that disrupts the slow col- 
lisional mode (operating in a Sweet-Parker-like geometry) (Cassak 
et al. 2005). This leads to a reconnection switch that triggers jet 
dissipation, which (as shown below) initiates prompt emission near 
the GRB jet photosphere. As mentioned in the introduction, the 
concept of a reconnection switch has been applied to many astro- 
physical phenomena, and in this paper we simply apply the same 
concept to GRB jets. 

The relevant length scale for fast collisionless reconnection in 
baryon-dominated plasmas is the proton skin depth 



d p = 



(1) 



for proton plasma frequency co pp = ^4m p e 2 /m p , proton number 
density n p , charge e, proton mass m p , and speed of light c. The 
relevant scale for pair-dominated plasmas is the pair skin depth 



(2) 



with pair plasma frequency u> pe = yjAnn eiol e 2 jm e , electron+pair 
number density n em = n e + n pairs , and electron mass m e . 

Collisions ensure that resistive MHD applies, which forces the 
current layer to avoid the fast Petschek-like geometry in favor of a 
Sweet-Parker-like geometry (see., e.g., Uzdensky & Kulsrud 2000). 
The Sweet-Parker solution requires pressure equilibrium across the 
current sheet. Assume, as valid for most of this paper, that the re- 
connecting field is not significantly weaker than the guide field. 
Then, the Sweet-Parker solution without a guide field can be used 
to estimate any quantities to order unity. The electromagnetic pres- 
sure (pem = "em, where uem is the electromagnetic energy density) 
and thermal gas pressure (which, as borne out in this paper, is dom- 
inated by photon pressure p y = u y /3, where u y is the photon energy 
density) balance via 



Pem ~ Py 



(3) 



Next, assume, as also borne out in this paper, that the collisional 
resistivity against current-carrying electrons and positrons is domi- 
nated by Compton scattering, then the corresponding magnetic dif- 
fusivity is 



n ~ (4/3)d 2 (u y cr T c)/(m e c 2 ), 



(4) 



(Goodman & Uzdensky 2008), for Thomson scattering cross sec- 
tion cr T . Then, for electromagnetically-dominated jets with rela- 
tivistic Alfven speeds, the collisional Sweet-Parker current sheet 
thickness is 



<5 S p ■ 



(5) 



where Lq is the comoving length of the current layer as estimated 
by the scale for variations in the electromagnetic field. 

The fast Petschek-like collisionless reconnection mode takes 
over when d p > <5 S p for baryonic-dominated plasmas, as discussed 
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above. Using the above equations, the condition for fast collision- 
less reconnection is then given by 

1< /^p_\ 2 (1 + n pii]l Jn e ) 2 
\Ssp) 8r y-sca /i 

where jl = iiEM/(PbC 2 ) is the electromagnetic energy per baryon 
rest-mass energy in the jet, p b » 2m p n p is the baryon rest-mass den- 
sity for a proton number density n p assuming protons and neutrons 
are at equal number densities, n e = n p is the baryonic-associated 
electron density from charge neutrality, T y sca as n e tot <r T L is the 
scattering optical depth, and L corresponds to the density scale- 
height for photons, electrons, and positrons for an emitting slab of 
length L (Popham & Narayan 1995). For a pair-dominated plasma 
in the limit that baryons play no role, the condition for fast recon- 
nection becomes 

de i - [ m ' Y i+n ^ ine) ( 7 ) 

~\Ssf) \m p ) 8T r-sca /i 

Future studies can clarify how fast reconnection operates in a 
plasma with baryons that typically dominate the rest-mass energy 
even if pairs have a non-negligible number density. 

Assume n pa j rs < n e , then the condition for transition to fast 
reconnection simply becomes 

r r , SC a <(8/2)" 1 . (8) 

Because MHD jets for GRBs have jl ~ 1 at large radii 
(Tchekhovskoy et al. 2010b), fast reconnection in GRB jets is pre- 
dicted to occur near the photosphere. Such a dissipation is required 
by dissipative photosphere models (see, e.g., Rees & Meszaros 
2005). Further, once fast dissipation starts, the value of jl drops, 
which forces the fast collisionless reconnection condition to be 
maintained. Even if fast reconnection were relativistic with speed 
order c, dissipation is still suspended until large radii so that a 
baryon-pure ultrarelativistic jet can form. 

The rest of this paper computes the radius where the recon- 
nection switch occurs, the photon optical depth at this radius to 
determine whether quasi-thermal photospheric emission is possi- 
ble, and some other observables. Realistic MHD GRB jet models 
are used since the polar jet's ideal MHD acceleration is substan- 
tially different than considered in prior reconnection jet models. 
The important pair contribution is non-trivial to compute because 
near the photosphere the photons and pairs are only marginally op- 
tically thick. Also, the plasma temperature near the transition to 
fast reconnection ends up low enough that pairs are deep in the 
suppressed regime. A simple exponential suppression factor would 
lead to more than an order of magnitude error in the transition ra- 
dius, so pairs must be treated more accurately. Despite pairs being 
in the suppressed regime, the pair number density can greatly ex- 
ceed the baryonic-associated electron number density and so cru- 
cially affect the transition radius as seen from Equation (6). We 
must also consider neutrino cooling near the compact object, be- 
cause strong cooling can lead to an effective drop in (5 SP (Uzdensky 
& McKinney 201 1) and so trigger fast reconnection causing the jet 
to dissipate before it is launched. Overall, this requires us to com- 
pute the properties of a jet filled with a complex of multiple recon- 
necting slabs, each consisting of baryons with electrons, pairs, and 
neutrinos at arbitrary optical depths. This is achieved by treating 
the current layer structure as similarly done for radiative accretion 
disks (see, e.g., Kohri et al. 2005). 




Figure 2. GRB jet structure near the central engine: A relativistic MHD jet 
is driven by a compact object with rotation frequency CI that causes the field 
line at foot point radius r[ p and opening half-angle 8f p to have a field line 
rotation frequency £2p (offset from the axis for clarity). At larger radii, each 
field line follows a collimating trajectory with opening half-angle 8{ , while 
the entire jet has an opening half-angle 6j corresponding to the largest 0{ p 
allowed by the central engine's accretion disk, corona, or wind. 

3 LARGE-SCALE JET STRUCTURE 

In this section, the large-scale axisymmetric structure of relativistic 
GRB jets from black holes or magnetars is presented. The large- 
scale jet structure acts as a key constraint on the reconnection 
physics by fixing the radial dependence of magnetic pressure. Dur- 
ing reconnection, the magnetic pressure balances the thermal pres- 
sure within current layers, which determines the contributions of 
baryons, photons, pairs, and neutrinos. 

Figure 2 shows the basic elements of the central engine and the 
production of an axisymmetric jet. The jet could be confined later- 
ally (up to a radius r mono ) by an accretion disk, disk corona, disk 
wind, stellar envelope, or some ambient medium. A generalized jet 
solution that applies for both small and large radii is obtained in 
Appendix A, which collects together results from recent analytical 
works and simulations (Tchekhovskoy et al. 2008, 2009, 2010b). 
However, the full solution can be cumbersome for obtaining sim- 
ple scaling laws for how the results depend upon model parameters. 

Therefore, this section outlines a simplified jet structure that 
approximates the full jet structure model and helps to highlight the 
basic elements of the full solution. The radius is assumed to be 
much larger than the size of the compact object (i.e. r » r fp ), 
at which a field line attaches at a footpoint radius (r fp ) and at a 
footpoint angle (9f p ) across the surface of the compact object (see 
Figure 2). The radius is also assumed to be larger than the decon- 
finement radius (i.e. r > r mmo ). The following expressions apply 
for < 9 < n/2 with an assumed equatorial symmetry. 

The degree of collimation inside r mono is given by a parameter 
v, such that for angles 9« 1 (i.e. cylindrical radius R <k r) one can 
show that the field lines obey 
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where the field collimates with 2 > v > out to r„ 



after which 



it follows the v = type radial geometry on a particular field line 
with fixed opening half-angle of Of giving 



ft* 



-v/2 



2 sin — 



(10) 



where the total jet opening angle (ft,) corresponds to the value of 
ft(ftp) for the largest value of ft p that is allowed by the presence of 
a confining medium near the compact object. 



For r > r mono » j-f p , the radial field strength is 



B r ~ #r,fp 



'IP 



(ID 



The ^-component is small and so can be neglected. The toroidal 
0-component is 

V V-1 / \-l 



2r fp f2 F Mr, 



tan(ft-/2), (12) 



where for rapidly rotating black holes or neutron stars rf p flp £ 
0.25c and Dp is the field rotation frequency one can set at each foot 
point. For a black hole angular rotation rate of SI = fi H = C/c)/(2r H ) 
with horizon radius r H an d dimensionless black hole spin j, the 
value of the field line angular rotation rate is Sip « Sl H /2. For 
Lorentz factor of y » 1, one can show that b 2 ss B^/y 2 , which al- 
lows one to obtain the electromagnetic pressure and energy density 
via pem = "em = b 2 /(8n), where \b\ is the comoving field strength. 
The baryonic rest-mass density is 



Pb ~ Pb,lp 



B, 

B r ,f P 



where y is some estimate of the Lorentz factor for r > r m 
define a magnetization parameter as 



B 



%np bSp c 2 



(13) 



We 



(14) 



which is used to define p b fp from S,. fp . Near the compact object, the 
value of £ is similar to the electromagnetic energy per particle. 

Overall, for a given 5,-.f p , r mono , ftp, and an estimate of y, one 
can determine the radial and angular dependence of b 2 , p b , Of, and 
6j, and other quantities defined in Appendix A; such as the electro- 
magnetic energy per unit rest-mass-energy (/}), the electromagnetic 
energy flux per unit rest-mass flux (p), the electromagnetic energy 
flux per unit mass-energy flux (cr), and the jet power (Pj). This 
MHD jet solution only has to be applicable up to the radius where 
significant dissipation occurs, because determining the dissipation 
radius is the primary goal of this work. 

For all models, we choose a rapidly rotating black hole foot 
point radius of r fp « 4.4km(M B H/M Q ) (M BH is black hole mass) 
with a rapid rotation rate such that r lp Sl F f v = 0.25c. The collapsar 
case with a black hole of mass M Bii = 3M Q is chosen. For the 
collapsar model r raono »3x 10 10 cm is typically chosen (i.e. radius 
of the progenitor star), while for the short-duration GRB compact 
object merger model r mono « 1.2 x 10 7 cm = 120km is chosen (i.e. 
extent of newly formed disk-corona-wind that helps collimate the 
jet ; Tchekhovskoy et al. 2008, 2010b). 

A black hole is chosen instead of a neutron star because it is 
more likely to be a generator of a powerful baryon-pure jet. First, 
the black hole cleans magnetic field lines of mass (MacDonald & 
Thorne 1982; Levinson & Eichler 1993), so the magnetization can 
be quite high in the black hole case. The magnetization p could 
be limited by neutron diffusion (Levinson & Eichler 2003), which 




Type A Type B & Type C Type D 

Figure 3. Generation of electromagnetic field reversals: Different jet field 
geometries at the jet base lead to current sheets with different orientations. 
Type A corresponds to a time-dependent polarity for an axisymmetric dipo- 
lar field. Type B corresponds to a time-independent non-axisymmetric mul- 
tipolar field. Type C corresponds to a time-independent axisymmetric mul- 
tipolar field. Near the jet base, types B and C are similar except the alter- 
nating field polarities are displaced either in the (^-direction (type B) or the 
^-direction (type C). Type D corresponds to a dipolar field that is unstable 
at large radii. All types are expected to some degree. 



still leads to quite high magnetizations of p ~ 10 3 — 10 4 (McKinney 
2005a) (corresponding to ( ~ 10 3 -10 4 ). Second, a black hole can 
produce a more powerful jet since its Kerr parameter is j ~ 1. For 
the magnetar case, the outflow only becomes highly magnetized at 
late time when the power has significantly diminished, and neutron 
stars have a Kerr parameter only up to j ~ 0.6 before break-up. 
These issues make it potentially more difficult for the magnetar to 
operate as both an efficient and powerful engine of a highly mag- 
netized ultrarelativistic jet (Metzger et al. 2007; Bucciantini et al. 
2008b; Metzger et al. 2010). 

One must use Appendix A to obtain an accurate dependence 
for all MHD quantities (i.e. including y) within the jet as a func- 
tion of radius and angle. Consider typical collapsar model param- 
eters that would lead to an ultrarelativistic jet with y ~ 1000. If 
<T = 10 4 , B„ fp = 3.2x 10 15 G, r mono = 3x 10 10 cm, and ft p = tt/2, then 
fi ~ 5400, y(r = 10 14 cm) « 800, Pj(r = r fp ) ~ 2.2 x 10 51 ergs -1 , 
[yff](r = 10 I4 cm) ss 18, and cr(r = 10 14 cm) a 6. At r = 10 14 cm 
there is about « 6 times less kinetic energy than electromagnetic 
energy, which can be tapped for prompt GRB emission if there ex- 
ists some mechanism to dissipate the energy. 



4 JET FIELD SUBSTRUCTURE: GENERALIZED 
STRIPED WIND 

In this section, we consider the process whereby small-scale field 
reversals and current sheets become embedded within the large- 
scale jet structure. The comoving length scale, Ao, of these jet field 
substructures plays a prominent role in later calculations because 
the dissipation rate due to collisional reconnection is dominated by 
the smallest value of Ao. 

Figure 3 shows the different types of field geometries that 
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are considered. All of these substructure types probably coexist 
to some degree. The jet structure solution presented in section 3 
was formally constructed with a single polarity for the electromag- 
netic field where the return field is presumed to connect back to 
the region slightly beyond the compact object. However, the polar- 
ity of the field can have arbitrary reversals without changing the 
large-scale force balance of the jet as long as the current sheets 
are inserted so that only the field direction changes. For an as- 
sumed field perturbation of approximately exp (i(k ■ Rj — £lpt)) with 
cylindrical radius Rj for the entire jet, Maxwell's equations give an 
overall smallest characteristic electromagnetic field length scale of 
A ~ l/(l/(7tRj) + n F /(27rc)). The comoving length scale (A ) gen- 
erally depends upon the substructure type. 

Substructure type A corresponds to a time-dependent alternat- 
ing polarity for a dipolar field near the compact object. The in- 
jection of field reversals occurs at an angular rate of ~ mflp near 
the compact object for some spherical harmonic quantum number 
m > 0. Such an injection can occur when the compact object's ro- 
tation axis is misaligned with the magnetic field axis as in the stan- 
dard striped wind (Coroniti 1990; Thompson 1994; Spruit et al. 
2001; Palenzuela et al. 2010), a neuron star's internal large-scale 
poloidal field dynamo flips the sign of the field, the compact object 
is fed alternating polarity from the accretion disk (Lovelace et al. 
1997), or a large-scale poloidal field accretion disk dynamo deter- 
mines the field near the black hole. Values of < m < 1 are allowed 
since a magnetic dynamo can operate on timescales much longer 
than the dynamical timescale. For example, an accretion disk could 
have m = 0.1 (Davis et al. 2010). 

Substructure type B corresponds to a time-independent non- 
axisymmetric (m > 1) multipolar field threading the compact object 
or accretion disk (Akiyama et al. 2003). The length scale of the 
electromagnetic field is determined by reversals at different foot 
point angles in the ^-direction near the compact object, and so type 
B is also related to the angular rate of ~ mQ. F . 

At large distances, substructure types A and B both have 
field reversals that are radially stacked such that the lab-frame 
A = 27ic/(m£lp). In the lab-frame, the pattern of reversals in A is 
advected to large radii at roughly the speed of light, so that the co- 
moving scale is A ~ Ay (Spruit et al. 2001). There may be many 
"stripes" separated radially by Ao. 

Substructure type C corresponds to a time-independent ax- 
isymmetric (e.g. spherical harmonic / = 1,2,3,... corresponding 
to dipolar, quadrupolar, etc.) multipolar field threading the com- 
pact object or disk. Notice that / cannot be zero because jet forma- 
tion requires a non-zero axial field at the jet base. The scale of the 
electromagnetic field is determined by reversals at different foot 
point angles in the ^-direction near the compact object. The field 
reversal scale is A ~ nRj/l. The scale of changes in the field is per- 
pendicular to the flow direction so that A = A ~ nRj/l. If yd > 1, 
then the jet becomes causally disconnected across in 9 due to finite 
wave speeds for any finite magnetization. The causally connected 
region has a size A ~ nRj/(ly9), so physical models effectively 
have / > yd when using the original A formula. The model de- 
scribed in figure 2 of Romanova & Lovelace (1992) corresponds 
to this substructure type C with 1 = 2 and with a prescribed struc- 
ture in (f>. Another example of type C is a magnetic tower model, 
where there is a single large current sheet that separates two re- 
gions (Lynden-Bell 1996; Uzdensky & MacFadyen 2006, 2007). 
This substructure is qualitatively different than substructure types 
A and B because Ao oc r while the other substructures do not de- 
pend upon r once y is roughly constant. 

Substructure type D corresponds to an ordered dipolar field 





Type A(detail) 



Type A and Type B 

Figure 4. Magnetic field reversals at large radii: Shows substructure types 
A-D at large radii where the toroidal field dominates and the flow carries 
many field reversals harboring current layers separated by length scale Ao. 
Type A corresponds to time-dependent reversals with Ao > y(2jrc/£lF). 
Type B corresponds to a "barber pole" of reversals with Ao ~ y(2nc / (mflp)) 
(m > is the toroidal spherical harmonic quantum number). Type C 
corresponds to reversals in 9 with Ao ~ (nRj/l) (I > is the poloidal 
spherical harmonic quantum number). Type D corresponds to reversals 
due to instabilities or turbulence at boundary layers and is some mixture 
of types A-C. In general, types A-D can be described by some particu- 
lar (more generally a power spectrum) of m > and / > to obtain 
A ~ WKicRj) + mSl F /{y2nc)). 

at the jet base, but the jet becomes unstable at large radii. Such 
instabilities could be due to intrinsic jet properties leading to the, 
e.g., current-driven magnetic kink instability or due to turbulence at 
boundary layers such as the jet- wind or jet-envelope boundary. For 
example, in the relativistically magnetized limit of MHD jets, if the 
comoving toroidal field exceeds the comoving poloidal field by fac- 
tors of a few to tens, then the growth rate for the kink instability is 
T ~ Rj/c generating structure from Ao ~ 2nR s to Ao ~ r (Narayan 
et al. 2009). Such kink modes would be stabilized by jet expansion 
or finite mass loading leading to causal disconnection across much 
of the jet except within 6 < 1/y (Spruit et al. 2001). In general, the 
scale of the perturbations could be randomly oriented across ~ 7rRj 
and generated on a time-scale of 2n/Q.p in the lab-frame. Hence, 
this substructure corresponds to some mixture of substructures A- 
C for m > and / > within A ~ l/(l/(nRj) + mQ. P /(y27Tc)). 

Figure 4 shows the jet substructures at large radii. Types A and 
B are similar, but in detail they differ. For types A-C, at all radii 
(including near the central compact object), the comoving poloidal 
field is at most order unity comparable in magnitude to the comov- 
ing toroidal field. For types A-C, the toroidal field dominates in the 
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lab and comoving frames for r » r mono . Generally, substructure D 
might correspond to having a strong poloidal field. However, mag- 
netic kink driven modes are expected to typically leave the jet's 
comoving toroidal and poloidal fields as comparable (i.e. when the 
so-called "safety factor" is of order unity in comoving frame), be- 
cause that is when the instability is initiated. Overall, one concludes 
that any guide field (which does not reconnect) in the current sheet 
is not much stronger than the reconnecting field. This means the 
analysis in section 2 and in later sections that assume a "weak" 
guide field are accurate to order unity or better, because the Sweet- 
Parker analysis is only significantly modified if the guide field is 
significantly stronger than the reconnecting field. 

In summary, the jet can develop various substructures that 
have comoving scale A satisfying 



1 / mflp 

A (/, m) nRj y2nc' 



(15) 



for some real numbers / > max(0, yd) and m > 0. Only dynamo 
and jet-medium interaction computations can determine (the power 
spectrum of) l,m, but many dynamos and instabilities saturate at 
low (/, m] ~ 1. The length scale (Ao) over which the comoving 
electromagnetic field changes is a natural scale over which any dis- 
sipation term enters into the comoving induction equation, and so 
A plays a central role in determining the comoving physics of re- 
connection discussed in the next section. At large radii, these modes 
correspond to the jet having current sheets with a guide field that 
is not significantly stronger than the reconnecting field, which al- 
lows one to estimate (to order unity) the properties of the collisional 
layer using the Sweet-Parker analysis without a guide field. 



5 RECONNECTION PHYSICS 

In this section, the conditions for a plasma to undergo collisional 
or collisionless reconnection in the comoving frame of the jet are 
computed. As described in section 4, the characteristic length scale 
of the electromagnetic field (before reconnection occurs) is given 
by A . As discussed below, such current distributions can be unsta- 
ble to collapse to a much thinner size d. The jet is assumed to be 
filled with a complex of multiple current layers each of thickness 

6 <k A . Instead of seeking a solution that is in global force balance 
within the jet, all fluid quantities (p, u, v h 5, in the background jet 
solution) are assumed to be roughly homogeneous within the cur- 
rent layer at any given radius. 

In section 5.1, the formation and length of current sheets is 
discussed, because the length enters into models of slow and fast 
reconnection. In section 5.2, the equation of state and emission 
rates for the current layer are summarized (see Appendix B for de- 
tails). In section 5.3, the transition from slow to fast reconnection 
is discussed (see Appendix C for more details). In section 5.4, how 
the electromagnetic field is modified by fast dissipation is com- 
puted. Finally, in section 5.5, the timescales for observed dissipa- 
tion/emission variability are estimated. 



5.1 Layer Formation and Length 

First, we must establish that current sheet formation is possible on 
timescales shorter than the jet propagation timescale to distances 
relevant for prompt GRB emission. Current layer formation starts 
with a distributed current of size A that spans a null point in the 
magnetic field strength. Such current distributions can collapse to 
current sheets on fast Alfven crossing time scales through a variety 



of mechanisms: 1) exponential collapse of the distributed current 
until the thermalized gas pressure within a thickness S < Ao can 
support the external electromagnetic pressure (Kulsrud 2005) ; 2) 
collapse as part of the evolution of the global magnetic field as oc- 
curs for the equatorial current sheet in a pulsar wind ; or 3) intrinsi- 
cally unstable collapse at an X-point (Dungey 1953, 1958; Imshen- 
nik & Syrovatskii 1967; Waelbroeck 1993; Loureiro et al. 2005). 
For example, X-type collapse forms current sheets on a timescale 



2v, 



■InS, 



(16) 



(Priest & Forbes 1986), where v A = c\b m / ^4n^\ ~ c is 
the Alfven speed outside the newly forming current layer, £ = 
p- m c 2 + u m + p m + b\j(4n), |£>j n | is the upstream comoving electro- 
magnetic field strength, w m is the upstream gas internal energy, and 
Pin is the upstream gas pressure. The relevant controlling parameter 
is the global Lundquist number, S = A v a /t], describing the ratio of 
the resistive time (r res = A 2 /rj for resistivity rj) to the typical elec- 
tromagnetic advection time (Ta = Ao/v A ). Typically S » 1 and 
so ideal MHD is valid on scales of order Ao. The above timescale 
shows that current sheet formation only depends logarithmically on 
the resistivity. For rapidly rotating black holes generating jets with 
7 ~ 100, half-opening angle 6j ~ 2°, substructure types A or B with 
m = 1, and S ~ 10 15 , current sheet formation occurs by r ~ 10' 3 cm 
as comparable with the expected radius of prompt GRB emission. 
Thus, current sheet formation is a relatively fast process that can 
occur during jet propagation. 

Second, the length of the current layer is needed because this 
is what enters the slow and fast reconnection models. The length 
of the current layer in the comoving frame (L ) can be estimated 
from the characteristic reversal length scale (A ) of the electromag- 
netic field, because this is the natural length scale that enters the 
MHD equations. The X-type collapse described above can feed off 
of perturbations introduced initially in the jet, producing structures 
with L a ~ A (Priest & Forbes 1986). In addition, the presence of a 
complex of multiple current layers can itself induce structures with 
L ~ A (Otto & Birk 1992; Yan et al. 1994). Ideal MHD type pro- 
cesses such as magnetic Rayleigh-Taylor instability can also pro- 
duce structure along the layer with L ~ A (Lyubarsky 2010a). 
Even in the absence of any current layers (i.e. m = and / = 1), the 
jet still has a fundamental characteristic scale of Ao(l, 0) = nRj. 

These arguments suggest that the length of the current layer 
should be within the range Ao < Lq < nRj. Because Ao ~ Rj is 
represented by I ~ 1 and for most of the paper A never appears 
explicitly independent from L , for simplicity, the allowed range 
for Lo is subsumed as an allowed variation in / and m. With the 
length of the layer determined, one determine the thickness of the 
slow collisional reconnection Sweet-Parker layer in section 5.3. 



5.2 Current Layer Pressure, Temperature, Density, and 
Emission Rates 

The reconnection physics is determined not only by the externally 
imposed magnetic field, rest-mass density, internal energy density, 
field reversal scale, and length of the current layer as described in 
the previous sections. In addition, the reconnection physics is deter- 
mined by the form of energy that the magnetic energy is dissipated 
into. 

The solution for the multi-species equation of state (EOS) and 
emission rates are described in Appendix B, which follows similar 
physics used for radiative accretion disks. In summary, the gas EOS 
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consisting of baryons, photons, non-pair-produced electrons, pair- 
produced electrons-positrons, and neutrinos is given by the total 
internal energy density 

U g = U b + U y + U e + Wpaijs + U v , (17) 

and total thermal gas pressure 

Pg = Pb + Py + Pe + Ppairs + Pv, (18) 

where each species is listed respectively and is a function of tem- 
perature (T) and baryon density (p b ). 

Given the proton number density (n p ), the baryon number den- 
sity (nt = pi, /m/,), charge neutrality such that the electron number 
density is n e = n p , then the electron fraction is Y e = n p /n b . A value 
of Y e ~ 1/2 is assumed (see Appendix B). Also, n e = n e - - n e +, 
where n e - and n e + are the total electron and total positron num- 
ber densities, respectively. Let the total electron+positron number 
density be n eM = n e - + n e + , then the pair-only number density is 

^pairs — Me,lot 

The reconnecting layer is treated as a radiating slab, with a 
total cooling rate (erg s -1 cm 4 ) for photons, pairs, and neutrinos of 

Qg = Qy + Gpata + Qv, (19) 

for each species respectively. The respective optical depths in the 
jet are given by T y , T pa j rs , and t v , which are separated into scattering 
and absorption optical depths (e.g., T r>sca and T rllbs , respectively). A 
non-negligible 2pairs occurs when T pa j rs < 1 as allowed when pairs 
leave the current layer by travelling along field lines. 

5.3 Transition from Slow Collisional to Fast Collisionless 
Reconnection 

In this section, we discuss how the current layers within the GRB 
jet make a transition from collisional to collisionless reconnection, 
which triggers the Petschek-like reconnection geometry leading to 
fast reconnection. Then, dissipation and emission from fast recon- 
nection can contribute significantly to prompt GRB emission. 

When the collisionless reconnection layer thickness is larger 
than the Sweet-Parker thickness (Sweet 1958; Parker 1957), the 
layer enters a Petschek-like reconnection geometry (Petschek 
1964) allowing for fast reconnection (see., e.g., Ma & Bhattachar- 
jee 1996; Uzdensky & Kulsrud 2000; Kulsrud 2001; Cassak et al. 
2005; Yamada et al. 2006). For current layers where plasma yS = 
Pg/Pb ~ 1, the condition that the plasma no longer obeys the resis- 
tive MHD equations on the scale of the current layer can be shown 
to be equivalent to 

6 S p> < dj = c/w pi , (20) 

where aj pi is the ion plasma frequency (Kulsrud 2005). Here, Ssv> = 
SspA ~' /2 , where Ssp ~ LqS~ >12 is the nominal Sweet-Parker thick- 
ness that is determined using a balance of gas pressure (p g ) and 
electromagnetic pressure (/?em) across the layer by following the 
Sweet-Parker analysis. The compression ratio A ~ max(l, Q g l 2sp) 
(where gsp is the dissipation rate in the Sweet-Parker layer) ac- 
counts for radiative cooling in the layer (Uzdensky & McKinney 
2011). Typically, A ~ 1 is found in this paper. The value of the 
Lundquist number, 5 , is obtained from Spitzer and Compton drag 
resistivities. For more details see Appendix C. 

When the above condition is satisfied, the reconnection oper- 
ates in a fast collisionless mode. The transition radius to fast colli- 
sionless reconnection is then given by 

rtrans = r[5 S P> = </(]. (21) 



This transition radius is used to identify where fast magnetic recon- 
nection starts and so leads to dissipation and emission. 

One interesting caveat is that once collisionless reconnection 
is initiated, the higher dissipation rate may lead to different plasma 
parameters such as the temperature. The loss of magnetic energy 
may raise the local temperature in the optically thick limit, but 
it can lower the temperature at slightly larger radii if the flow is 
already optically thin because a loss of magnetic energy leads to 
lower gas pressures. If resistivity increases due to dissipation, then 
the Sweet-Parker thickness increases and may lead to recovery of 
the collisional mode of reconnection. In such a situation, the plasma 
may enter into a "marginally collisionless" reconnection mode in 
which <5sp' = dj is maintained as the plasma dissipates and con- 
tinues to enter and exit the fast and slow reconnection modes. It is 
difficult to estimate the average rate of reconnection in such a sce- 
nario (Uzdensky 2007b; Goodman & Uzdensky 2008). However, in 
some studies, the transition back to slow reconnection was shown 
to be avoided by a hysteresis effect (Cassak et al. 2005). 

Slow collisional reconnection can be as slow as Sweet-Parker, 
giving an inflow velocity of v r ~ v A S~ l/2 <K v A ~ c, or it may be 
as fast as v, ~ 0.01i>a (Uzdensky et al. 2010). Various studies have 
found that the fast collisionless reconnection rate is of order v r ~ 
0Av A or even v r ~ v A (for details see Appendix C4). For simplicity 
and definiteness, the reconnection rate is chosen to be v r ~ 0.1c in 
all cases. The exact number used for v r only moves the expected 
dissipation radius and does not significantly affect the transition 
radius or any other calculations. Once significant dissipation has 
occurred, the value of v A drops below c, but in all cases we will find 
it only drops to v A ~ 0.01c by the time the slow collisional mode 
would recover. Also, since v A dropping below c indicates most of 
the electromagnetic energy flux has been dissipated (i.e. the goal of 
this work), we simply avoid detailed discussions about the solution 
at large radii where this occurs. 



5.4 Effect of Magnetic Reconnection on the Jet Structure 

In this section, we determine how the jet electromagnetic energy 
density, b 2 /(87r), vs. r is modified by fast dissipation once the tran- 
sition to fast reconnection has been triggered. The resulting func- 
tion for b 2 vs. radius modifies the original large-scale jet structure 
dependence given in section 3 and Appendix A. The radial extent 
over which significant dissipation occurs also determines the pho- 
ton/pair/neutrino opacity integrals used to determine the equation 
of state and emission rates for the radiative current layer in Ap- 
pendix B3. 

The electromagnetic energy dissipated in the jet is given by 
the electromagnetic enthalpy flux into both sides of a current layer 
where each side harbors a region of size A /2 from which magnetic 
flux can be accumulated for each sheet. The comoving electromag- 
netic dissipation energy density rate is then 

d !b 2 \ b 2 v r 
dt m \ An A 

where v r is the comoving rate of reconnection, corresponding to 
the inflow of magnetic flux into the current layer. Notice that the 
relevant length scale is Ao, the distance between layers instead of 
the length of the layer, because magnetic flux must be brought in 
from the space between layers. 

For the jet structure, the relativistic MHD energy-momentum 
equations of motion are dominated by the radial electromagnetic 
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field advection term' for either an electromagnetically-dominated 
jet with 7 » 1 or for y and Of roughly constant as occurs for r > 
'"mono (as applicable to all astrophysically relevant models studied in 
this paper). Then, Equation (22) is a sink term to the quasi-steady 
relativistic MHD energy equation given by 



= - 



1 



r 2 sin 



d r (r 2 sm6b 2 y 2 c)-4nQ EM y. 



(23) 



For Qem — > (i.e. v r — > 0) corresponding to the ideal MHD jet 
solution, the comoving electromagnetic energy density is a power- 
law dependence with b 2 oc r~ 2 . In general, Equation (23) can be 
solved for b 2 (r) to give 

b 2 (r) 



b 2 (r ) 



■■ exp 



dr 



2v r 



cyA Q 



(24) 



which assumes ty oc v A ~ c as relevant for electromagnetically- 
dominated jets. Once va <k c, significant dissipation (i.e. the goal 
of this work) has already occurred, and then we simply avoid de- 
tailed discussions about the overall solution for such radii rather 
than seeking the correct b 2 vs. r in the weakly magnetized regime. 

For m + modes, the value of A is fixed (i.e. does not scale 
with r for r > r mono ), and then the jet solution has an exponential 
term such that, 



b 2 [r = r^] 



exp - (r - r ttans ) 



2v r 
cyA 



(25) 



for a radius of transition to fast reconnection r ttans and r > r trans . 

On the other hand, for I > modes A oc r, and then the 
dissipation does not lead to exponential decay and instead leads to 
a modified power-law given by 



b 2 [r — Ttrans] \r tri 



(26) 



for r > r trans . In general, one obtains r^ ss /r ttsas ~ f-nwOj/Vivr) f or 
a fraction / of magnetic energy dissipated (e.g. r diss ~ 10V trans for 
/ = 5, v r /c = 0.1, y6j = 10, and / = 1/2). A simple conclu- 
sion would be that negligible change in b 2 occurs for large enough 
x r = (y6j/l)(7Tc/v r ) as due to a given /-mode becoming causally 
disconnected across the jet. However, the / modes generate plas- 
moids with radial size Lo that may remain fixed in size once sig- 
nificant dissipation occurs. This alternative assumes, perhaps opti- 
mistically, that for / modes the magnetic flux can be accumulated 
and dissipated from all directions instead of only in the 6 direction. 
So instead of assuming some I modes fail to dissipate, an alterna- 
tive is also considered that the / modes produce plasmoids similar 
to substructures A and B except with an effective Ao — > Lo whose 
value is fixed with radius as evaluated at the transition radius. 

Once fast reconnection starts at r aws , then Equation (25) 



This calculation can be compared to that in Lyubarsky & Kirk (2001); 
Kirk & Skjaeraasen (2003b), and also compared with Drenkhahn (2002) 
who made some simplifications to obtain their equation 31 giving their 
equation 37. This simplification was also used to obtain equation 1 in 
Drenkhahn & Spruit (2002b). Drenkhahn et al. concluded that no more than 
50% of the electromagnetic energy flux can be converted to optically thin 
radiation and that the rest goes into kinetic energy flux due to jet acceler- 
ation. However, no works have yet fully included radiative cooling that (in 
the optically thin regime) can remove thermal energy on a timescale com- 
parable to the timescale for heating and jet acceleration. For simplicity, we 
assume dissipation only leads to radiation (as might occur in a fully radia- 
tive calculation) but does not lead to jet acceleration. This issue does not 
affect any of our interpretations of results. 



shows that reconnection completes a single e-folding of dissipa- 
tion/emission by 

ycA 



'"dis; 



'"trans 



2v r 



(27) 



assuming m modes or the optimistic case for / modes. For example, 
consider a model with m = OA, y ~ 300, a rapidly rotating BH 
such that Q. F ~ 0.25c/r H , A as 2ncy/(mQ. F ) as 3y x 10 8 cm, and 
v r ~ 0.1c. This model gives r diss as r ttans + 10 14 cm, which results 
in a reasonable radius for prompt GRB emission for small enough 
'"trans- For r > r diss , most of the dissipation has completed and leaves 
relatively few photons or pairs. So, this radial extent can be used 
to determine the opacity integral for photons, pairs, and neutrinos 
within a relativistic jet as done in Appendix B3. 

Between the transition and dissipation radii most of the jet's 
electromagnetic power, P {l r M) , is lost as radiation from plasma 
blobs moving at Lorentz factor y being observed at some angle 
away from the jet axis. For a lab-frame radial dissipation range of 
dr ~ ycA /v r , lab-frame cross-sectional area ~ ttR 2 , and lab-frame 
dissipation rate as yQ EM , the total isotropic jet luminosity is 

Lj ~ (yQ EM )(ycA /v r )(irR 2 ) ~ y 2 b 2 cR 2 ~ pf M) as P p (28) 

which shows that the lab-frame loss rate of electromagnetic energy 
equals the lab-frame electromagnetic jet power. This also shows 
that the reconnection rate, v r , does not determine the luminosity 
from the region dr. Instead, a larger number of layers within the 
range dr dissipate simultaneously for smaller v r /c, so a fixed jet 
power leads to a fixed total isotropic luminosity. 



5.5 Variability from Fast Reconnecting Layers 

In this section, we compute various timescales from reconnecting 
layers that may be imprinted in observations of the prompt GRB 
emission. The jet substructure type determines the temporal struc- 
ture of the prompt GRB emission via the number of dissipating 
layers. Also, emission timescales can be linked to the total event 
timescale, the event timescale per unit number of reconnecting cur- 
rent layers embedded in the jet, the dissipation timescale of a single 
current layer, and the thermal photon radiation emission timescale. 

Consider a reconnecting layer (plasma blob) that begins to 
dissipate at a transition radius r tmns with a comoving dissipation 
timescale of dt m . This is related to the lab-frame timescale via 
dhd, as ydt co . Photons emitted at the start of fast reconnection travel 
outward as the plasma blob also travels outward. The observed pho- 
ton time difference from beginning to end of current sheet dissipa- 
tion is dt bs as dtub/(2y 2 ), such that the observed photon timescale 
accounting for light-travel effects is dt obs as dt co /(2y). 

For the m modes and for the optimistic case for I modes, the 
comoving dissipation timescale is dt iiss _ co ~ A /(4v r ), correspond- 
ing to an observed timescale of 



dt a \ 



Ao 

Syv r ' 



(29) 



Each layer present in the jet over the duration of the event 
gives one expected emission timescale. For substructure types A 
and B the lab-frame cycle time to produce a single slab is dt\-^ m ~ 
2n/(mQ.p). Therefore, one expects a variability timescale of 

2n 



dt \ 



mQ. F ' 



(30) 



associated with each shell that emerges from r = r trans . For sub- 
structure type C, there are I simultaneous emitting layers that in 
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principle could dissipate stochastically in time rather than going off 
simultaneously. In that case, for an event duration T, one expects a 
variability timescale of 

T 

dt obs j ~ —. (31) 

Reconnection is expected to proceed stochastically, but it will 
also have a short timescale associated with the transit time of fluid 
through a single reconnection layer of exhaust length Lq. The co- 
moving transit time is order rff co ,transit = Lo/va, such that each re- 
connection event has an observed time scale of 

L 



dt 0] 



2yv A 



The thermal photon radiation emission timescale is 



(32) 



(33) 



corresponding to the timescale over which the thermal radiation 
contributes fully to its component of force balance against the mag- 
netic pressure in the current layer complex. 

For scale-height H and mean free path of A = \/(dr/ds), the 
photon diffusion timescale is 

"fobs.d 



2r 



H 



Ac 



(34) 



where for H = H a t, s one has H/A = t, the optical depth. The ob- 
server absorption diffusion timescale is denoted dt ohwi(i when set- 
ting H = H abs , while the observer scattering diffusion timescale is 
denoted rff bs,sd when setting H = // sca . 



6 RESULTS: JET DISSIPATION 

This section computes several results obtained with the procedure: 

• Given the independent variables (, B r [p , v, r mono , and 9 {p and 
the jet structure from section 3 (specifically, using Appendix A), 
one obtains p b (r, 9), b 2 (r, 6), y(r, 0), Q. ¥ (r, 6), and 6 f {r). 

• Given Z, and m and the jet substructure from section 4, one 
obtains the single current layer spacing of A (/, m) and layer length 
of Lo(l, m) at any point in the jet. 

• The multiple current layer complex equations determined in 
section 5, with the non-ideal MHD jet solution for b 2 (r) (equa- 
tions 25 and 26) used after r = r tmm , are iteratively solved as a 
set of simultaneous equations for the temperature T, pair number 
density n pa ir S , and the compression ratio A. A single (unique) solu- 
tion is always found. 

• Given T, n pairs , and A, remaining quantities are iteratively 
computed. 

The solution at each radius is obtained independently 2 . The transit 
time through the current layer is roughly the equilibrium timescale, 
which is sufficiently small compared to the jet flow time for this 
to be a valid approximation except for radii much larger than the 
transition or dissipation radii. 

In what follows, a fiducial model is discussed in detail in sec- 
tion 6.1, a model parameter exploration is presented in section 6.2, 
the simplified jet structure given in section 3 is used to write down 
scaling laws for results as functions of independent variables, a 
comparison is made to prior work in section 6.4, and finally in 



Each solution takes 1-24 hours to compute on a modern Xeon core, 
where over all models about 100, 000 CPU-hours are required. 



section 6.5, a discussion is provided about how the reconnection 
switch model applies to other (non-GRB) jet systems in order to 
identify whether GRB jets are special. 



6.1 Results: Fiducial Model with Full Jet Structure 

Consider a fiducial model of a rapidly rotating black hole with a 
magnetization f = 10 4 that leads to 100 < y < 1000 ; a radial mag- 
netic field at the compact object of B nfp = 3.2 x 10 15 G that leads to 
Pj ~ 10 50 -10 52 ergs -1 ; a collimating field geometry with v = 3/4 
as roughly consistent with a shocked stellar envelope (see discus- 
sion in Tchekhovskoy et al. 2008, 2010b) ; a deconfinement radius 
of r mono = 3 x 10 10 cm as roughly the radius of a Wolf-Rayet star 
in the presupernova phase ; and a foot point opening half-angle of 
f?fp ~ n/2 as applicable to a neutrino-dominated accretion disk that 
takes on a quite thin geometry near a rotating black hole where the 
jet is launched (see, e.g., Kohri et al. 2005; Chen & Beloborodov 
2007 and for simulations of such magnetized thin disks see, e.g., 
Shafee et al. 2008; Penna et al. 2010). The chosen f = 10 4 corre- 
sponds to an electromagnetic energy flux per unit rest-mass flux of 
p ss 5400 and an electromagnetic energy per baryon rest energy at 
the field line foot point of p = b 2 /(8np b c 2 ) ss 17000. For simplic- 
ity r fp fi Fi fp = 0.25c is assumed, as applicable for a rapidly rotating 
neutron star or black hole. 

Given the ambiguity in whether small / (with in = 0) modes 
actually efficiently dissipate (see end of section 5.4), let the fidu- 
cial model correspond to / = 1 (dipolar jet with no separated cur- 
rent sheets in the ^-direction), and let m = 0.1 (low order non- 
axisymmetric mode). The choice of m = 0.1 is motivated by recent 
studies of accretion disks showing non-axisymmetric modes gener- 
ated on tens of dynamical times (Davis et al. 2010). 

This fixes the values of all the free model parameters, defining 
our fiducial model explored in this section. The figures show 40 
positions in radius along a field line within the jet. The effects of 
varying the model parameters on the solution are explored in the 
next section. 

Figure 5 shows the structure of the jet along a single field line 
from the compact object to large radii. At a radius of r ~ r mono , 
the opening angle rapidly becomes constant and the Lorentz factor 
jumps up due to a rarefaction wave that forces the jet towards the 
monopolar solution. The jet is cold in the sense that pem <; P s 
and y » u g /(p b c 2 ) up to the dissipation radius, so the cold ideal 
MHD jet structure solution is roughly valid at large radii even when 
current sheets are dissipating. 

For MHD jets, the Lorentz factor is limited such that strictly 
y < p at all radii. The Lorentz factor is dominated by radial motion 
for our jet solutions. As discussed later, the loss of electromagnetic 
energy occurs near the jet photosphere, so no dominate thermal fire- 
ball is generated. This loss of electromagnetic energy causes the to- 
tal energy flux per unit rest-mass flux (p) to drop down to the bulk 
Lorentz factor. 

Ideal MHD jets are characterized by starting with a poloidal 
field stronger than the toroidal field at the foot point on the com- 
pact object, while at larger radii the toroidal field dominates. In 
the present model, the field near the compact object is in the 
super-critical QED regime. The comoving field at large radii is 
\b\ ~ \B$\ly, \b\ » \B r \ before significant dissipation occurs because 
B r oc r~ 2 in the monopole regime beyond r = r mom , and y remains 
roughly constant. In the ^-direction the electromagnetic pressure is 
constant as required by force balance. 

Figure 5 also shows how the ideal MHD jet solution is mod- 
ified by the current sheet dissipation. Beyond the transition radius 
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Figure 5. Jet structure for fiducial rapidly rotating black hole model with f = 10 4 , B r f v = 3.2 X 10 15 G, v = 3/4, r mono = 3 X 10 10 cm, 0f p = n/2, and modes 
m = 0.1, / = 1, for quantities along the outer-most field line. Top-left panel: Opening half-angle {Of) as thick solid line. The location of r = r mono is shown by 
the thin vertical line. Top-right panel: Lorentz factor (y) as thick solid line. Total energy flux per unit rest-mass flux (// > y) as mostly horizontal dotted line. 
In this and other panels, the label "Id" marks the same quantity for that line type, but corresponding to the ideal MHD jet solution without electromagnetic 
dissipation. The ideal MHD jet solution has constant fi and y continues to rise beyond the dissipation radius. Bottom-left panel: Magnetic field components 
(B r as solid line and -Bp as dotted line). Ideal MHD jet solution at large radii follows roughly -B^ oc 1/r. Bottom-right panel: Baryon rest-mass energy 
density (pbC 2 ) as solid line and comoving electromagnetic energy density (fo 2 /(8/r)) as short-dashed line. The thin horizontal solid line shows \b\ = ^qed, the 
QED critical field. As with the toroidal field strength, the comoving electromagnetic energy is dissipated away in the non-ideal MHD jet solution beyond the 
transition radius. This is in contrast to the ideal MHD jet solution that would continue with roughly b 2 oc r~ 2 , which is a similar power-law scaling as rest-mass 
density at large radii. Overall, by r a 10 14 cm the jet obtains an opening half-angle of 1°, a Lorentz factor of y as 800 with y0j a 18, and a jet power of 
Pj ; as 8x 10 51 ergs" 1 . 



at r aim s ~3x 10' 3 cm, dissipation proceeds due to fast collisionless 
reconnection. This causes \b\ (and so |S^|) to drop exponentially for 
several e-foldings. Fast dissipation finally ceases once the transi- 
tion back to slow collisional reconnection occurs due to low tem- 
peratures leading to a dominant Spitzer resistivity. Despite the sig- 
nificant drop in the comoving toroidal field, only by r ~ 6x 10 15 cm 
does it become smaller than the comoving poloidal field. Beyond 
this radius, any current sheet analysis requires consideration of a 
guide field, but this is beyond the dissipation radius so that all of 
our work remains valid. 

Figure 6 shows the temperature, collisional and collisionless 
current layer thicknesses, number densities, and optical depths for 



neutrinos, photons, and pairs as functions of r. The sequence of 
events from small to large radii is as follows: 1) The optically 
thick current layer's temperature decreases until pairs drop out 
r « 10 10 cm ; 2) Photons dominate the gas pressure. But after the 
photon absorption photosphere at r « 2x 10 n cm, the photon densi- 
ties decrease faster ; 3) The loss of photons causes the temperature 
to begin to increase at r ~ 2 x 10 12 cm, which eventually leads 
to a reemergence of pairs at r « 4 x 10 12 cm (with n pllirs = n e at 
r x 8 x 10 12 cm) and then to a slight quenching of the rise in tem- 
perature because the pairs suspend the decrease in the photon scat- 
tering opacity ; 4) The transition to fast collisionless reconnection 
occurs at r ttms »3x 10 13 cm, where n pairs ss ln e ; 5) The transition 
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Figure 6. Reconnection layer solution using the jet model shown in Figure 5. Top-left panel: Current layer temperature (T). The vertical line shows the 
location of r mono , where the jet becomes deconfined. The jumps in temperature are resolved by several points and explained in the text. Top-right panel: 
Fast collisionless (Petschek-like) ion skin depth (dj) as solid line and slow collisional (Sweet-Parker-like) layer thickness (iSspO as short-dashed line. Notice 
that dj = Ssf at r a, 3 X 10 13 cm, where fast collisionless reconnection is triggered. As in Figure 5, the label "Id" marks the case when the ideal MHD jet 
solution is used to determine the properties of the collisional layer. Bottom-left panel: Electron number density (n e ) as solid line, pair number density (n pa i rs ) 
as short-dashed line, and neutrino number density (n y ) as long-dashed line. Jet expansion leads to cooling and a loss of pairs at r ~ 10 10 cm until they reemerge 
at r as 8 X 10 12 cm. Bottom-right panel: Neutrino scattering optical depth (t„_ sc11 ) as solid line, photon scattering optical depth (r r-sca ) as 

1 at 



with « pairs 

short-dashed line, and pair scattering optical depth (r pi 
r as 5 X 10 13 cm. Overall, the key result is that the fast collisionless reconnection mode is triggered at r ■■ 
coincides with where the jet becomes optically thin for photons. 



a ) as long-dashed line. The thin horizontal line shows an optical depth of unity, where r r 



3 X 10 13 cm where ; 



1, so the transition roughly 



radius is nearly coincident with where photons begin to free-stream 
leading to T 7iSca ~ 1 at r ss 5 x 10 13 ; 6) Eventually the pairs drop 
out again atr« 10 15 cm due to the pair absorption and scattering 
optical depths approaching unity, after which the temperature rises 
rapidly again ; and 7) An e-folding of electromagnetic dissipation 
occurs by rji ss ss 2 X 10 l3 cm, leading to a decrease in temperature 
due to a loss of electromagnetic and hence thermal pressures. 

The key point from Figure 6 is that d, = <5 SP < at /- trans ss 
3x 10 13 cm, marking the transition radius to fast collisionless recon- 
nection. Beyond r trans , the collisional resistivity is no longer effec- 
tive and dissipation proceeds due to collisionless processes that op- 
erate on the larger scale of the ion skin depth. This disrupts the cur- 
rent layer into a Petschek-like geometry, which allows a relatively 



fast reconnection rate of ty ~ 0.1-1.0c rather than a slower Sweet- 
Parker-like rate. This is why we identify this as a catastrophic re- 
connection switch. 

Note that the radial range where electrons are degenerate is 
far inside the transition radius, so degeneracy does not affect the 
resistivity nor the value of r trans . Also, note that even with the loss 
of pairs for r ~ 10 10 cm-4 X 10 I2 cm, the proton-electron plasma 
remains highly collisional, which forces the transition to fast colli- 
sionless reconnection to be at large radii. 

The transition to fast collisionless reconnection might lead to 
a marginally collisionless state that forces dj ~ <5 S p< . Using the full 
non-ideal MHD jet solution shows that the collisional layer thick- 
ness rises once dissipation occurs, and fast dissipation ceases when 
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the two thicknesses become roughly equal because the collisional 
mode recovers. However, by this radius, significant electromagnetic 
dissipation has already occurred - and in fact this is what allows the 
collisional mode to recover. In addition, the possibility for a return 
to collisional reconnection was tested by using different v r up to 
v r ~ c, which leads to a different radial dependence of b 2 . We find 
that marginal collisionality was unable to be established for any v r 
until there has already been significant electromagnetic dissipation. 
Also, notice that the dissipation included in the model does not it- 
self control the transition to fast reconnection. This was tested by 
using the ideal MHD jet solution for b 2 vs. r, which shows that the 
thickness for the collisional layer remains much smaller than for 
the collisionless layer once fast collisionless reconnection would 
have started. 

The condition for fast reconnection in electron-positron plas- 
mas is the pair-dominated switch condition d e = <5 S p< . This condi- 
tion occurs at r ~ 6 x 10 15 cm in the ideal MHD jet solution that 
assumes no dissipation up to this radius. If the baryon-dominated 
switch condition applies and dissipation did occur, then the pair- 
dominated condition is reached at r ~ 2 x 10 15 cm after signifi- 
cant dissipation has already occurred. However, the pair-dominated 
switch condition is not expected to be relevant for the fiducial 
model or for a broad range of variations in model parameters. Pairs 
only reemerge after the photon scattering opacity becomes of or- 
der unity, and the transition to fast collisionless reconnection has 
occurred to order unity. Also, the total mass-energy of the plasma 
is found to be carried by protons by the radius where d e = (5 S p' or 
di = S S p' ■ So, the condition for fast reconnection may still be dom- 
inated by the protons, although future studies should consider how 
fast reconnection proceeds when pairs dominate in number density 
but ions dominate in energy density. 

Figure 6 also shows the scattering optical depth for neutri- 
nos, photons, and pairs. The neutrinos are never optically thick in 
this model. The quantity t 7 corresponds to the opacity towards the 
observer accounting for locally generated pairs and downstream 
baryonic electrons as computed via Equation (B22). The photons 
become optically thin to scattering at r ~ r lrans , which is consis- 
tent with our basic argument presented in section 2. This means 
that the reconnection switch mechanism occurs once non-thermal 
and quasi-thermal photons can be produced. Beyond r ~ r trans , the 
non-ideal MHD jet solution has a photon scattering opacity that 
deviates only slightly from that of the ideal MHD jet solution. So 
dissipation is not crucial for lowering the photon scattering opti- 
cal depth at larger radii. The pair scattering opacity is dominated 
by Coulomb interactions when the temperature is low, but then the 
pair density is negligible (e.g. beyond r ~ r iiss leading to a jump in 
the pair opacity). For the ideal MHD jet solution, the pair optical 
depth roughly follows that of the photons and continues to drop at 
all larger radii. 

Perhaps the most distinct non-intuitive feature of the solution 
is the two-step rise and final drop in temperature. Let us consider 
a simplified calculation to explain why the temperature is forced 
to rise once the jet becomes optically thin to scattering. Cyclo- 
synchrotron radiation dominates over free-free emission at all radii, 
so the optically thin energy density loss rate is 



Q ~ n em -a T c(y e - l)w EM , 



(35) 



where y e as l + & e /(T - l) (recall & e = k B T/(m e c 2 )) is the electron 
thermal Lorentz factor. Then, T r-al)s as gAo/(cw , r ) is the absorption 
optical depth from Kirchhoff 's law, and the scattering optical depth 
is T r sca ~ <x r rc pairs A since pairs (when present) typically dominate 



the photon scattering opacity. The energy density of photons is ob- 
tained from the two-stream approximation using Equation (B14), 
which in the limit r r , a bs ^ 1 and r r , sca ~ 1 gives 



u y « 3.2T ria b s «o,i 



(36) 



where u 0y = a m(i T 4 is the optically thick radiation energy density, 
and p y = Uy/3 is the radiation pressure that dominates over other 
pressures. Then, the pressure equilibrium condition across the cur- 
rent layer of p EM = w EM = p y determines the temperature implicitly 



2.9 



r- 1 



+ 1.4 



0, 



r- l 



(37) 



when T 7 , a b S < 1 (otherwise the temperature is set by p , r = Pem)- 
So a drop in the scattering optical depth leads to a rise in the tem- 
perature in this limit. One also obtains 



"EM 



(38) 



which indicates that the absorption opacity drops significantly 
when the electromagnetic energy density is not high enough to 
produce optically thick photons. Equation (37) shows that once 
the flow becomes optically thin, the temperature is forced to rise 
up to T ~ 10 9 K regardless of the temperature in the optically 
thick regime. The full solution has a slightly lower temperature of 
4 x 10 8 K because the cyclo-synchrotron emission is slightly domi- 
nated by higher energy electrons than those with the mean temper- 
ature. This gives a solution closer to using T = 4/3 than T = 5/3 
in the above simplified equations. Overall, the behavior of the tem- 
perature vs. radius can be explained as follows: 1) The temperature 
decreases with radius due to adiabatic expansion in the optically 
thick regime ; 2) The temperature rises when the photons become 
optically thin due to the need to maintain pressure equilibrium with 
electromagnetic field ; 3) The slight flattening in the temperature 
during its rise (two-step rise) is because pairs reemerge and increase 
the scattering opacity, which suppresses the rapid temperature rise ; 
4) Eventually the pair opacity also decreases and the pairs drop out, 
which leads to a continued rise in the temperature as the scattering 
opacity decreases ; and 5) Finally, the temperature decreases inde- 
pendently of the photon scattering opacity once a significant loss of 
electromagnetic energy has occurred. 

The scattering and absorption photospheres lead to thermal 
emission and might correspond to the observed peak energy. In ad- 
dition, low energy electrons' cyclo-synchrotron emission is self- 
absorbed and then Comptonization can dominate their emission 
leading to a Comptonized thermal component as the peak energy 
corresponding to Epeak ~ 10keV-2MeV (Rees & Meszaros 2005; 
Thompson et al. 2007). Electromagnetic dissipation at the tran- 
sition radius (coincident with the photon scattering photosphere) 
leads to a thermal peak of £ pcak sss yT ~ 20MeV for the fidu- 
cial model value of y ss 800. The thermal peak obtained at the 
photon absorption photosphere is smaller at T ~ 10 7 K, giving 
£ P eak ~ 600keV at this relatively high y ss 800. This peak energy is 
roughly consistent with the observed peak energy for cosmological 
long-duration GRBs. 

Figure 7 shows the Lundquist number, observer timescales, 
energy density loss rates, and pressures within the collisional cur- 
rent layer at each radius. 

For the ideal MHD jet solution, one would estimate that the re- 
sistivity is dominated by Compton drag at all radii, with the Spitzer 
(due to proton-electron Coulomb collisions) becoming nearly as 
important at r ~ 7 x 10 8 cm and at r ~ 2 x 10 15 cm. For the full 
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Figure 7. Continued reconnection layer solution (otherwise identical to Figure 6). Top-left panel: Lundquist number (5) as short-dashed line. As in Figure 5, 
the label "Id" marks the ideal MHD jet solution. Using the ideal MHDjet solution for the current layer shows that the Lundquist number stays roughly constant, 
while the full non-ideal MHD jet solution shows the Lundquist number drops as associated with a high Spitzer resistivity that dominates over the Compton drag 
resistivity. Top-right panel: Observer time for transit through the current layer (A bs,tra) as solid line, observer time for single e-folding for fast reconnection 
dissipation (dr b s ,diss) as short-dashed line, observer time for m = 0.1 radial structure induced variability (dt ob!iJn ) as dotted line, and observer time for photon 
emission from the current layer (dt obSiy ) as long-dashed line. Overall, the observer timescales during significant electromagnetic dissipation are order Is and 
span 0.001s to 10s. Bottom-left panel: Photon energy density cooling rate (Q y ) as solid line, neutrino energy density cooling rate (Q v ) as long-dashed line, and 
upper limit for cooling rate in strong-cooling regime (2sp) as dotted line. The layer is not in the strong-cooling regime because Q s < Qsp, which means the 
layer thickness can be estimated by the standard Sweet-Parker type solution. Bottom-right panel: Photon pressure (p y ) as solid line, pair pressure (p pa irs) as 
short-dashed line, and neutrino pressure (p v ) as long-dashed line. Baryonic-associated pressures are always negligible until significant dissipation has already 
occurred. The photon pressure is comparable or dominates all other sources of thermal pressure. 



dissipative MHD jet solution, Spitzer resistivity begins to play a 
role only after significant electromagnetic dissipation occurs once 
there is a significant decrease in the temperature within the layer at 
r ~ 6x 10 15 cm. The associated Lundquist number is order S ~ 10 20 
and eventually drops at large radii due to MHD dissipation. 

Figure 7 also shows various observed timescales for photon 
emission. Once fast collisionless reconnection kicks in at r ttal i S , 
these timescales give an estimate for observed variability. The rise 
in the timescales at large radii is due to the jet being optically thin 
and so the diffusion timescale reaches the light crossing time that 
increases with radius. The transit time for fluid to pass through 
the current layer would lead to variable emission directly due to 



reconnection, while some of the other timescales give variabil- 
ity due to photon diffusion and emission rates. The m mode pro- 
duces many pulses due to new reconnection layers passing be- 
yond rtrans. Emission would decay over the timescale for com- 
pleting dissipation. The observer time for photon absorption diffu- 
sion from the current layer photosphere (dr b s , a d) and the observer 
time for photon scattering diffusion from the current layer pho- 
tosphere (dr obs . s d) roughly follow df obs , r . Between r ttans and r dlss , 
these timescales are df obs ,tra ~ 0.002s-0.007s, dt obsm ~ 0.004s, 
*>bs,ad ~ 0.001s-0.04s, dt obsM ~ 0.002s-0.04s, dt ob%y ~ 0.004s- 
0.06s, and rff bs,diss ~ 0.005s. After a few e-foldings of dissipation, 
the timescales are within the range of ls-lOs. These timescales are 
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roughly consistent the observed GRB pulse durations that typically 
are ~ 0.5s, but range from 0.01s to 10s (Norris et al. 1996). Be- 
cause the timescale is always associated with a fixed radius, pulse 
timescales do not evolve over the event duration as happens in the 
internal shock model. 

Figure 7 shows that neutrinos dominate the energy density 
loss rate up to r ~ 10 10 cm, after which photons dominate. Near 
the jet base, neutrino energy density loss rates do not exceed gsp, 
which means that the current layer does not compress due to radia- 
tive cooling and the non-radiative Sweet-Parker analysis is approx- 
imately valid. The gas energy cooling rate Q g ~ <2sp within the 
dissipation region after fast collisionless reconnection is triggered, 
so radiative cooling has become marginally dynamically important 
to the current layer structure. 

Figure 7 also shows that the photon pressure is comparable or 
dominates other pressures at all radii including in the optically thin 
regime. Electron and baryon pressures are negligible at all radii. 

Lastly, consider some checks on the composition of the jet. 
The electron and positrons are found to be non-degenerate at 
all radii except for a region between r mono and about 100r mono . 
Nucleons are computed to everywhere be non-degenerate and B- 
equilibrium implies Y e = 1/2 everywhere as consistent with this 
work's assumptions. However, the nuclear statistical equilibrium 
(NSE) timescale can be compared to the jet flow time such that 
complete NSE occurs at cnse ~ cydtxsE- For the fiducial model, 
the value of r NS E is always found to be much greater than the local 
radius. NSE only completes to about a tenth of a single e-folding 
NSE timescale. So, the composition of the jet would primarily be 
determined by the accretion disk (either directly by advection and 
reconnection or by neutron diffusion) rather than by NSE or B- 
equilibrium. 

6.2 Results: Parameter Dependence for Full Jet Structure 

In this section, several model parameters are explored to determine 
how they affect the results. All fiducial parameters are held fixed 
except a single parameter that is allowed to vary. For each model 
parameter and across all radii, a total of 40x40 solutions (i.e. radius 
x model parameter) are sought. 

Figure 8 shows results for B r (r fp ) vs. r at f = 10 4 for which 
ji ss 5400. We find that a stronger electromagnetic field near the 
compact object leads to a transition to fast reconnection at larger 
radii. For expected collapsar parameters with B r{p ~ 10 15 G, the 
transition to fast collisionless reconnection occurs at r ~ 10 13 - 
10 14 cm. In all cases, at r > 10 14 cm the Lorentz factor is y ~ 800 
for the ideal MHD jet solution. The electromagnetic energy flux per 
unit mass-energy flux is cr ~ 6, indicating a large reservoir of en- 
ergy exists to be electromagnetically dissipated. This solution has 
y6j ~ 18 for all field strengths at a radius of r ~ 10 14 cm, while 
asymptotically yOj ~ 20. Jet breaks can occur because y6j > 1. 
For all of parameter space, A ~ 1 so the current layer is not 
in the strong-cooling limit. The break in the behavior of r diss at 
B r ,fp ~ 10 14 G is due to the Lorentz factor saturating at large radii. 

When t„ > 1 , nuclear statistical equilibrium is found to hold 
and free nucleons are generated in the jet independently of the 
accretion disk composition. Above B r j p ~ 10 15 G, near the com- 
pact object nuclear statistical equilibrium partially holds and the 
jet contains a non-negligible fraction of free nucleons that typically 
freeze-out before they might convert to a-particles at large radii. 
Below S rfp ~ 10 13 G, the jet is found to never be in nuclear statisti- 
cal equilibrium on jet flow timescales, and the composition would 
be set by the accretion disk that feeds the jet with baryons. The 
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Figure 8. Contour plot for radial magnetic field strength at the field line foot 
point (B,-(rf p )) vs. r showing the location of transition radius for fast colli- 
sionless reconnection (r trims as left-most diagonal solid line bounding cross- 
hatched region), location of a single e-folding of fast collisionless dissipa- 
tion (rdi ss as right-most diagonal solid line bounding cross-hatched region), 
the location of the photosphere defined by the optical depth to radial infinity 
T y = 1 (diagonal dotted line roughly following the transition radius), re- 
gions where electrons dominate pairs in number density (region with shade 
lines going from bottom-left to top-right corner), the region where neutrinos 
are optically thick (region in upper-left corner), and the region where the co- 
moving electromagnetic field is QED super-critical (left of the short-dashed 
line). For this and similar figures, note that the pairs drop out again within 
the dissipation region due to significant dissipation and drop in temperature. 
However, for clarity this is not shown here because it always starts within 
the cross-hatched region. Overall, the jet becomes optically thin towards 
the observer near the transition radius to fast collisionless reconnection and 
typically long before dissipation completes, which allows significant dissi- 
pation to produce non-thermal or quasi-thermal emission. 



/3-equilibrium value of Y e ~ 1/2 holds for all field strengths. How- 
ever, NSE timescales are typically long compared to the jet flow 
time, so reaction rates need not be considered. Thus, for the field 
strengths required to explain typical cosmological GRBs, the com- 
position in the accretion disk would determine the jet composition 
rather than equilibrium conditions or reaction rates. 

Plots similar to Figure 8, but with different f, show similar 
results. For all f > 10 2 , the flow is optically thin when fast col- 
lisionless reconnection initiates. Consider the case with f = 10 2 
for which fi = 55. The transition radius to fast collisionless recon- 
nection is also roughly linear in log-log and varies from r ams = 



10 57 cm 
B, 



r H to r tt 



10 cm for a variation of B r fv = 10 G to 



r,fp = 10 G, respectively. Also, r iiss ~ 10r, rans for B r[p < 10 G, 
while otherwise r^ ss ~ r ttans . In addition, the region where pairs 
drop out shifts down along with this r ttans line. Asymptotically, 
7 ~ 50 and cr ~ 1. However, the solution always has y8j < 2, which 
implies jet breaks may not be discernible for such solutions even if 
they produce efficient non-thermal or quasi-thermal emission. For 
collapsar parameters and this value of ( ~ 10 2 , the transition to fast 
collisionless reconnection occurs at r ~ 10 15 cm. 

Figure 9 shows how changing the effective magnetization pa- 
rameter, f, determines the solution for otherwise typical collapsar 
parameters that we adopted for our fiducial model in section 6.1. 
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Figure 9. Contour plot for the magnetization (f) vs. r showing ydj = 1 
(lower mostly horizontal dashed line) above which jet breaks can occur, and 
otherwise like in Figure 8. Neutrinos are never optically thick for this case. 
By r > 10 14 cm and for f ~ 10 4 , the Lorentz factor is y ~ ,u/7 ~ ( 1 13. Also, 
cr ~ 1 for ( = 10 2 and increases for increasing f, while y0 ~ 1 for ( = 10 2 
and increases for increasing f . This implies jet breaks are predicted for a 
broad range in f . For f > 5000, the flow becomes optically thin towards the 
observer near the transition radius to fast collisionless reconnection. This 
shows for collapsar parameters with B r ~ 10 15 G that the transition to fast 
collisionless reconnection occurs at r ~ 10 14 cm iff ~ Iff corresponding 
to n ~ 5400 and y ~ 800. 



This shows that dissipation starts at r ~ 10 14 cm if £ ~ 10 4 , which 
gives y ~ 800 with cr ~ 6 and y0 ~ 18. We consider this as a fidu- 
cial model that is fairly insensitive to other free parameters. This 
type of progenitor system generates a jet fast enough to avoid the 
compactness problem, harbors about 6 times more energy available 
for prompt emission than for afterglow emission (assuming the en- 
ergy is mostly dissipated and radiated away instead of accelerating 
the jet), and can produce jet breaks. For £ < 5000, the dissipa- 
tion completes within the region where photons scattering opacity 
is unity, which forces a non-negligible fraction of photons to ther- 
malize. Such events should have a thermal emission component in 
their spectra. In this parameter space for the ideal MHD jet solu- 
tion, the comoving toroidal field has increased to a few times the 
comoving poloidal field, triggering the kink instability (Narayan 
et al. 2009), at r ~ 10 133 cm for £ > 10 35 and at r ~ 10 108 cm 
for f = 50. However, for all f, the transition to vigorous kink in- 
stability still occurs beyond r mono where y6j » 1, so dissipation 
takes place only near the rotational axis within 6 < l/y due larger 
angles being causally disconnected. This angle is typically small 
compared to the jet opening angle, and so this region contains very 
little electromagnetic energy flux due to both the electromagnetic 
power dependence on angle and the efficient ideal MHD accelera- 
tion at small angles (Tchekhovskoy et al. 2009, 2010b). 

All other free model parameters were similarly explored for 
otherwise fixed fiducial parameters. The free parameter v controls 
the field geometry ranging from monopolar (v = 0) to parabolic 
(v = 1) with the most likely value for the collapsar model being 
v ~ 3/4. The solution dependence on v = ( 10~ 4 , 2) is weak, where 



r trans varies from 10 13 cm to 10 14 cm for v = 10 4 to v = 2 and where 
r diss varies from 10 131 cm to 10 15 cm for v = 10~ 4 to v = 2. 

The free parameter r mono controls when the flow transi- 
tions from a collimating solution (e.g. as caused by a stellar en- 
velope or accretion disk corona/wind) to a monopolar solution 
as occurs when the jet becomes unconfined (i.e. no collimating 
agent is present). We find that for r mono = {GM/c 2 , 10 18 cmj, 
the transition radius to fast collisionless reconnection increases 
slowly/monotonically from (r trans , r mono ) of (10 13 , 10 57 )cm to 
{10 15 , 10 18 )cm. For r mono > 10 15 3 cm, jet breaks are not expected 
because yOj < 1 once r ~ 10 15 4 cm. 

Consider short-duration GRBs formed by BH-NS (black hole 
- neutron star) or NS-NS mergers that involve only an accretion 
disk instead of both an accretion disk and a stellar envelope. Then, 
''mono ~ 10 7 cm is applicable, because this corresponds to the cylin- 
drical radial extent of the accretion disk that collimates the jet. 
Such a system is found to have a jet opening half-angle of order 
8j ~ 0.4 x 23°, and the transition radius to fast collisionless recon- 
nection is at r trans « 10 13 cm. 

The dependence on 6 fp = (0. 16, tt/2), over which much of 
the jet power resides, is very weak. The transition radius to fast 
collisionless reconnection only varies from (10 13/7 cm to (10 138 cm. 
There is little electromagnetic energy flux for field lines with # fp < 
0.1, so dissipation there is not important compared to the total jet 
dissipation. 

Figure 10 and Figure 1 1 show that the variations in / and m 
lead to significant changes in the behavior of the solution. A sig- 
nificant portion of the dissipated electromagnetic energy is ther- 
malized and produces a hot magnetic fireball when / > 10 2 and 
m > 0.8 rather than direct non-thermal emission. This effect is 
caused by the larger /, m leading to a smaller current sheet length 
and hence a smaller Sweet-Parker scale 5 S p, which then becomes 
equal to the ion skin depth at smaller radii. One may generally 
have a spectrum of /, m modes, and then dissipated electromagnetic 
energy contributes both thermal and non-thermal components to 
the photon spectrum. A stronger constraint on the I modes is also 
shown using Equation (26) with / = 0.5 corresponding to dissipat- 
ing half of the electromagnetic energy flux. This would be the dis- 
sipation radius if the collisional structures (that formed before or at 
the transition radius) continue to expand with radius as naively pre- 
dicted by A . It seems more likely that once the collisionless mode 
is triggered, then the dissipation decouples from any collisionally- 
induced structures, in which case the weaker constraint can be used. 

Overall, for collapsar type parameters, the parameter space 
study shows that the transition to fast reconnection often occurs 
beyond or near the photosphere unless the magnetization is too 
low (f < 10 3 ). The transition occurs near, but below, the photo- 
sphere when the / mode is too high (/ > 10 2 ) or the m mode is too 
high (m > 1). For collapsar type parameters, the transition radius 
is roughly linear in a log-log plot of the foot point field strength 
(B,-(r fp )[G]) vs. radius, such that, by coincidence, the transition ra- 
dius happens to have a value in centimeters roughly the same as the 
value of the foot point field in Gaussian units. Interestingly, the nat- 
ural m ~ 1 spiral mode (McKinney & Blandford 2009) demarcates 
a boundary where the transition to fast reconnection occurs at the 
jet photosphere. So small parameter variations can lead to events 
with thermal or non-thermal components in the prompt GRB spec- 
trum. 
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Figure 10. Contour plot for the substructure / mode vs. r and otherwise like 
in Figure 8. Unlike small /, large / causes the transition radius to fast col- 
lisionless reconnection to occur within the optically thick region. The line 
labeled by r dissalt corresponds to Equation (26) for / = 0.5 corresponding to 
dissipating half of the electromagnetic energy flux. This is a stronger con- 
straint on the dissipation radius, which shows that for I < 100 the dissipation 
is very slow even for fast reconnection due to the (possible) continued radial 
and lateral expansion of the jet on dissipating structures. 
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Figure 11. Contour plot for the substructure m mode vs. r and otherwise 
like in Figure 8. Unlike small m, large m > 0.8 causes the transition radius 
to fast collisionless reconnection to occur within the optically thick region. 

6.3 Results: Parameter Dependence for Simplified Jet 
Structure and Current Layer Physics 

In this section, the simplified jet structure presented in section 3 
is used to obtain formulae that describe the essence of the above 
results. First, one must obtain an estimate of the bulk jet Lorentz 
factor, y, such as can be obtained from the full jet structure so- 



lution in Appendix A. Second, one must approximate the current 
layer physics, which is possible for a some portions of the relevant 
parameter space. The results are not generally accurate to within an 
order of magnitude, but the calculation at least highlights the basic 
ideas and principles. 

For typical collapsar parameters and at the transition radius 
to fast collisionless reconnection, the photon pressure dominates, 
and the absorption opacity (dominated by cyclo- synchrotron pro- 
cess) is very small while the scattering opacity is not small. Cyclo- 
synchrotron radiation can be approximated with non-relativistic 
electrons having a delta function distribution at energy 3k]$T '/2. Ra- 
diative cooling within the layer is often weak, so that A = 1 can be 
chosen. The Coulomb logarithm is order 30. Photon drag resistivity 
generally dominates over the Spitzer resistivity up to and far be- 
yond the transition radius. A non-relativistic temperature approxi- 
mation is accurate, although pairs do provide non-negligible cyclo- 
synchrotron radiation and opacity. The pair and photon opacities 
are such that g[T pairs ] and g[r y ] (see Appendix B14) are roughly 
comparable. These approximations greatly simplify the force bal- 
ance condition given by Equation (C7), leaving a couple of equa- 
tions for T and rc pairs to be solved numerically. 

Even with these simplifications, a difficulty in obtaining a 
closed-form approximation is that the temperature suppression of 
pairs is transcendental. However, one can leave n pairs as a free pa- 
rameter, which can be estimated from the full solution or the ap- 
proximate solution in the preceding paragraph. 

We find that the value of f is implied from the value of y, 
because at large radii y ss f /14 for models that end up with <x order 
several. For models near the collapsar type parameters and f < 10 6 , 
the expressions in the rest of this section tend to be accurate to a 
couple orders of magnitude. For very large ( > 10 6 , the dependence 
of y on free parameters can become difficult to estimate without the 
full jet structure, and so the below expressions should be used with 
caution when arbitrarily varying y. 

For compactness of the below formulae, we rescale the free 
model parameters by collapsar values: f fp = r fp /r g with r g = 
3GM e /c 2 , r mono = r mono /(3 x 10 10 cm), B r , fp = B r , fp /(3.2 x 10 15 G), 
f = r /700, fp = sin [f9 fp /2]/ sin [tt/4], and % = <T/10 4 . Both / and m 
modes are considered but in general have different properties. For 
the case considered, however, n pa j rs = « pa j rs /(6n c ) is reasonable for 
both mode types. Such a choice for n pairs is a reasonable estimate 
for only a limited range of model parameters such as chosen for the 
fiducial model near the transition radius as estimated by this sim- 
plified model. Generally, both y and n pairs need to be chosen based 
upon the results from the full solution. Otherwise, for / modes, free 
model parameters are rescaled as follows: r ; = r/(l x 10 14 cm), 
and L j = (Lo)/(7ri?j et /Z) such that I = y6j, where L 0J assumes that 
typically L scales with /fj et as when m = 0. For m modes, free 
model parameters are rescaled as follows: r m = r/(l x 10 14 cm), 
and L Q m = (L )/(yc/(mQ. F )) with m = 0.1, which assumes that L 
does not scale with R sa as for I = 1 and m t 0. As applicable to 
most of parameter space, we assume the optical depth integrals use 
L ~ A < rjy for m modes and L ~ A > rjy for / modes. 

Then, the transition radius is 

rtrans,i ~ W cm U r{f L QJ r fp r mono y fp y n pairs C, , (;Sy) 
for I modes, and 

W ~ 10 14 cm B rSp I£ ^ y- 112 T" 4 , (40) 

for m modes. This transition radius is consistent with our full solu- 
tion from the previous sections. These estimates are valid for both 
non-relativistic or relativistic temperatures, whereas other scaling 
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laws obtained later tend to require a non-relativistic approximation 
for reliable scalings. Notice that the transition radius is only weakly 
sensitive to n pa i rs . In the case where ji pairs < n e , one can roughly 
take rc pairs — > n e in the above expressions. For remaining estimates, 
"pairs » n e within the current layer is assumed, although one can 
readily solve the system of equations that emerge in the other limit. 
The dissipation radius is 

fdissj ~ fW., + 10 14 cm (B, fp L ,,) 4/3 r^ 2 "' 3 f 2 ^ 9 2 £ (n palrs l)^ f~ 

(41) 

for / modes, and 

F diss.m ~ f trans,m 

+ 10 14 cm L , m f fp f, (42) 

for m modes. 

The optical depth towards the observer for thermal photons is 

T y , ~ 2 U,i h H P ai rs J X + 0.04 B% p f 2 p v r mmo f- 2 fj l t\ (43) 

for / modes, and 

T y , m ~ 1 L , m ffp ;j palrs + 0.1 B 2 fp f] v v F mmo y- 1 f' 1 l~\ (44) 

for m modes. The left term corresponds to the optical depth due 
to the current layer, and the right term corresponds to the optical 
depth due to electrons associated with baryons downstream in the 
jet. The local pair contribution tends to dominate the downstream 
electron contribution to r y . However, as L decreases, the right term 
dominates. 

Optically thin synchrotron could be the source of the observed 
high-energy photons near the peak of the energy spectrum for 
prompt GRB emission. Internal shocks are a fast heating process 
that could accelerate particles to high non-thermal energies. How- 
ever, reconnection is a slow heating process such that heating and 
optically thin synchrotron cooling must be in balance, which leads 
to low-energy particles and too low synchrotron photon energies to 
explain the observed peak energies (unless, e.g., high-energy elec- 
trons cool downstream where field is much weaker) (Ghisellini & 
Celotti 1999). 

Quasi-thermal emission could instead explain the observed 
prompt spectrum (Rees & Meszaros 2005; Thompson et al. 2007; 
Giannios 2008). There is necessarily incomplete photon thermal- 
ization because r y < 1 for r > r trans . Equation (37) already es- 
tablished that at a unity scattering optical depth there is a stable 
temperature of T « 4x 10 8 K independent of model parameters. Be- 
cause the transition to fast reconnection occurs once T y sca = 1, this 
provides a naturally stable temperature. The temperature is high 
and would only explain the observed peak energy for y < 50. How- 
ever, several scatters may be required to thermalize the photons, 
and incomplete thermalization might lead to lower temperatures 
allowing up to y ~ 1000 while still obtaining consistency with the 
observed peak energy (Thompson 1994; Giannios & Spruit 2005; 
Giannios 2006; Giannios & Spruit 2007; Giannios 2008). 

Incomplete thermalization at the scattering photosphere may 
lead to a temperature closer to that at the absorption photosphere. 
Let ® e = ® f /0.0004, and assume n pa j rs <k n e as applicable for many 
model parameters, then the temperature at r y _ abs = 1 satisfies 

©^Pp) 2 + O.O 5 (^p) 6„ (45) 

for both / and m modes, such that for fiducial parameters ® e = 1 
or & e = 0.0004. Typically, the left-hand term dominates, then & e ss 
(^/(B rfp f fp )) 419 . Then the observed thermal peak energy is at 

£ P eak,obs ~ 7 ® c 140keV, (46) 



which is close to the observed peak energy. However, a more de- 
tailed study is required to determine whether £ pea k oc y 13 ' 9 is too 
much variation to be consistent with, e.g., the Ghirlanda relation 
(Ghirlanda et al. 2004). 

One of the shortest timescales for observed prompt variability 
can be estimated from the transit time (as seen by an observer) for 
plasma to pass through the reconnection layer as given by 

*obs,tra,i ~ 0.02 s L 0J n f- 2 , (47) 

for / modes, and 

0.004 s L„,,„ f^, (48) 

for m modes. This and the other timescales in section 5.5 give a 
range of values from 0.001s up to about 10s for both m and I modes, 
which is roughly consistent with observations. For field substruc- 
tures types A, B, and C, this suggests that the prompt GRB emis- 
sion probes the turbulent field generated in the accretion disk near 
the black hole. For field substructure type D, the prompt emission 
probes the turbulent layer that develops between the jet and stel- 
lar envelope or accretion disk corona/wind. While one expects re- 
connection to proceed stochastically, the saturated dissipation of a 
collection of current sheets would lead to a delay in the dissipa- 
tion of current sheets just passing the fast transition radius. So one 
expects a correlation between pulse widths and intervals, and one 
expects no spreading of pulses at later times because the transition 
radius is roughly fixed. More modeling is required to validate these 
suggestive agreements with observations (Piran 2004). 

6.4 Comparison with Other Works 

Pioneering studies of reconnection in GRB jets by Spruit et al. 
(2001) determined the radius (their equation 61) at which dissipa- 
tion proceeds. In this present paper's language, they assumed sub- 
structure type A (and type B with m = 1) to obtain a dissipation 
radius of r d ss (nc) / (efl^y 2 , where e ~ 0.1 controls the reconnec- 
tion velocity normalized to Alfven velocity. They assumed that fast 
reconnection occurs over the entire length scale Ao, giving a lab- 
frame timescale for dissipation of T ~ y(A /v r ). They conclude 
that r d ~ 2 x 10 12 cm for £1 F w 10V and a fixed value of y w 300. 

One should check that this condition holds at smaller radii 
where y ~ 1 in order to ensure the jet can form in the first place 
without dissipation disrupting the jet formation process. At small 
radii, the jet has y w 1 + R£l F /c and 9 w (r/r fp r v/2 (2 sin (0 fp /2)). 
For any v, one can solve for the dissipation radius accounting for 
the changes in the Lorentz factor. Their conclusions appear accurate 
as long as the jet contains substructure type A (or B with m ~ 1) 
and dissipation occurs over A . Some problems may arise in their 
model because the dissipation radius decreases inversely with m. 

More serious is that if the jet contains field reversals in 9 as in 
substructure type C, then their equation (57) shows that dissipation 
would occur near the jet base because y9j ~ 1 at small radii for 
a magnetized jet 3 . GRMHD simulations of accreting black holes 
show that substructure type C, which is generated by MHD turbu- 
lence, is quite common. This suggests that only by having a slow 
reconnection rate at small radii and a fast reconnection rate at large 
radii can one ensure magnetized jet formation in the first place. 
The reconnection switch model presented in this paper naturally 
provides such a mechanism. 

3 A hot MHD jet can be pushed to slightly larger ydj at small radii. 
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6.5 Applications to Other Jet Systems 

AGN and x-ray binaries may also exhibit a transition to fast col- 
lisionless reconnection. Here we consider only a couple example 
systems. The ranges r g < r mom < 10 1( V S and 1 < £ < 10 6 are 
investigated. One expects ( ~ y < 100 (Fender 2003). 

First, consider the case of the jet in M87 (similar arguments 
apply to Blazars), which has a black hole with mass M « 6x 10 9 M o 
accreting at M ~ 10~ 2 M o /yr with a bolometric luminosity of 
L bo] ~ 10 42 ergs~' (Gebhardt & Thomas 2009), and so operates 
at a radiative efficiency of order r\ < 0.01 corresponding to a low- 
luminosity radiatively inefficient accretion flow. Equipartition argu- 
ments combined with GRMHD simulations of accretion disks and 
the Blandford-Znajek funnel region (McKinney 2004) then imply 
that the magnetic field strength at the base of the jet is of order 
B r ,f P ~ 10 3 G. 

Our solutions show that the M87 jet is already in the fast colli- 
sionless regime at the base of the jet (see also Giannios et al. 2010). 
This suggests that an electromagnetically-dominated jet from M87 
would not survive unless an ordered dipolar field is present at the 
jet base (McKinney & Blandford 2009). Because the magnetic field 
should be well-ordered near the black hole to produce a jet in the 
first place, this suggests that pairs may dominate the mass-energy 
because baryons would only enter the jet base via reconnection or 
turbulent diffusion. This may have important implications for polar- 
ized disk and jet emission and composition (see., e.g., Dexter et al. 
2010; Broderick & McKinney 2010; Shcherbakov et al. 2010). 

Second, consider the case of the jet in the BH x-ray binary 
GRS1915+105, which has a black hole with mass M ~ 14M Q 
(Greiner et al. 2001) accreting at near Eddington rates during out- 
burst and production of a fast transient jet (Fender et al. 1999). The 
system has M ~ 2.4 x 10~ 8 M Q /yr (Kording et al. 2006) and a bolo- 
metric luminosity of L bo i ~ 10 38 ergs~', so operates at a radiative 
efficiency of order r\ ~ 0. 1 corresponding to a radiatively efficient 
accretion flow. Equipartition arguments combined with GRMHD 
simulations of accretion disks and the Blandford-Znajek funnel re- 
gion (McKinney 2004) then imply that the magnetic field at the 
base of the polar jet is B r fp ~ 10 9 G. 

Interestingly, for GRS 1915 + 105 with £ < 10 4 and only within 
20 gravitational radii, the jet is in the slow reconnection regime (see 
also Goodman & Uzdensky 2008). Such a confined region where 
reconnection is essentially marginally collisionless could promote 
the production of the fast transient jet as triggered by a transition 
from slow to fast collisionless reconnection. More modelling is re- 
quired to investigate this effect. 



7 DISCUSSION 

For GRB jets, why must reconnection be slow enough near the cen- 
tral engine and yet fast at large radii? 

Fast relativistic reconnection at speeds up to order v r ~ c 
might occur (Lyutikov & Uzdensky 2003; Lyubarsky 2005). If 
this operated near the central engine and the jet contained current 
sheets, then dissipation of the jet would occur in situ and no jet 
might emerge or would be heavily baryon-loaded (Beckwith et al. 
2008; McKinney & Blandford 2009). If somehow a rate of v r ~ c 
was suspended until y » 1, at which point relativistic time dila- 
tion would stall dissipation, then still the dissipation radius would 
be inside the photon's photosphere (i.e. r diss < r trans , where T y ~ 1 
at r = ftjans). For example, for m = 0. 1 and otherwise fiducial pa- 
rameters, having v r ~ c leads to r iiss ~ 10 13 cm while T y = 1 at 



r = 10 14 cm. All dissipation would complete inside the photosphere 
and non-thermal emission would be unlikely. Hence, in general, 
fast reconnection should be delayed somehow. 

On the other hand, slow collisional reconnection leads to a de- 
layed dissipation rate relative to fast reconnection by a large factor. 
For example, in the most extreme case of Sweet-Parker reconnec- 
tion, the slowdown is by a factor of ~ Vs" , where S » 1 is the 
Lundquist number (see Equation 29). At the transition radius, fast 
reconnection has an observed timescale on the order of a second 
for the fiducial model. A typical Lundquist number at the transition 
radius is 5 ~ 10 20 . This indicates that a time delay of about 10 10 s 
(i.e. hundreds of years) would occur before the slow reconnection 
can dissipate a sizable fraction of the available energy, which is 
an unacceptable delay. Thus, slow collisional Sweet-Parker recon- 
nection (with any jet substructure type) cannot be responsible for 
powering the prompt GRB emission. Therefore, fast reconnection 
is required, but under general considerations it must be delayed un- 
til large radii to avoid complete thermalization of the photons. The 
reconnection switch mechanism satisfies both these requirements. 

An alternative to fast collisionless reconnection is fast colli- 
sional reconnection. Collisional reconnection may proceed faster 
than Sweet-Parker due to secondary tearing instability (plasmoid 
dominated reconnection) or MHD turbulence. Preliminary simu- 
lations show that, for otherwise fixed parameters, non-relativistic 
Petschek may only be 10 times faster than non-relativistic colli- 
sional reconnection (Loureiro et al. 2007; Kowal et al. 2009a; Sam- 
taney et al. 2009; Cassak et al. 2009; Loureiro et al. 2009; Huang 

6 Bhattacharjee 2010; Uzdensky et al. 2010). While fast collision- 
less reconnection would dominate this collisional mechanism once 
the plasma becomes collisionless, one should consider the effect at 
smaller radii. The jet starts with velocity v ~ c on scales of order 
the black hole horizon and the Lorentz factor grows rapidly with 

7 oc Rj. Unless the reconnection rate were v r — > c, the relativistic 
time dilation stalls reconnection until y flattens-out for r > r mono . 
Assuming fast collisional (e.g. plasmoid-dominated) reconnection 
is 10 times slower than the typical speeds of fast collisionless re- 
connection, then the dissipation radius is at 10 times the distance 
given by Equation (41) and Equation (42). So, fast collisionless 
dissipation already completes before fast collisional reconnection 
even becomes important. Also, for the fast collisional plasmoid- 
dominated mode, it remains unclear whether causal disconnection 
across the jet would allow for plasmoid chains to grow. As with 
suppression of the magnetic kink instability, turbulence may be 
avoided except within the narrow region of 6 < 1/y where the jet 
is still causally connected. Because this region contains very little 
power, this would lead to negligible dissipation compared to the 
total jet dissipation. 

Interestingly, GRB systems are quite unique compared to 
AGN and x-ray binaries. GRBs exhibit a transition from slow colli- 
sional to fast collisionless reconnection at roughly billions of grav- 
itational radii from the central compact object. On the other hand, 
X-ray binary systems may exhibit a transition to fast reconnection 
within tens of gravitational radii, while AGN jets tend to always be 
in the fast collisionless regime. This means that while GRB systems 
might sustain a disordered field at the jet base, it appears unlikely 
that a disordered field can be sustained at the jet base in most AGN 
and some x-ray binary systems. Significant dissipation and emis- 
sion can then occur when the jet harbors current sheets. This sug- 
gests that ultrarelativistic jets may be more difficult to obtain from 
AGN and x-ray binaries unless the magnetic field geometry con- 
sists of an organized dipolar field (McKinney & Blandford 2009). 
For AGN, such a requirement is tolerable if the observed coherence 
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length of magnetic fields in the ISM, which is high enough to trap 
a significant flux near a supermassive BH, is assumed to be typical 
of galaxies with AGN (Narayan et al. 2003). 

This suggests that current sheets within jets in AGN and x- 
ray binary systems may limit how efficiently electromagnetic en- 
ergy flux can be converted into kinetic energy flux through ideal 
MHD acceleration. Normally the electromagnetic energy flux per 
unit mass flux (41, measuring the degree of mass-loading) is under- 
stood to limit the terminal Lorentz factor. However, dissipation of 
electromagnetic energy in current sheets can lead to significant de- 
creases in fi. This may explain why GRB jets tend to have higher 
Lorentz factors than jets from AGN and x-ray binary systems. 

Another interesting result is that the current-driven kink mode 
is an unlikely source of current sheets. As discussed in section 4, 
the \m\ = 1 kink instability operates when the magnitude of the co- 
moving toroidal field (|Zty|) is about 3-10 times the magnitude of 
the comoving poloidal field (\b p \) (Narayan et al. 2009). The co- 
moving toroidal field's relative growth only occurs for r » r mono . 
In the fiducial model for the ideal MHD jet solution with 6 fp = tt/2, 
these comoving fields become comparable by r ~ 10 13 cm and the 
comoving toroidal field is about 10 times larger by r ~ 8 x 10 13 cm. 
Across all models, \b$\ only grows to 3- 10 times \b p \ for r » r mono . 
So, only after this radius would one expect vigorous comoving kink 
instabilities. However, beyond r = r mono , the jet loses significant 
causal contact across in angle due to the typical value of y6j ~ 10- 
20. So, one only expects the region within 9 < 1/y around the 
rotational axis to undergo the instability. For typical parameters of 
yOj ~ 15, this region contains negligible electromagnetic power. 
Further, the structure of the jet has \b$\ <K \b p \ at such small angles 
due to the even more efficient acceleration that occurs closer to the 
rotational axis (Tchekhovskoy et al. 2009, 2010b). Overall, kink in- 
stabilities appear ineffective at leading to significant dissipation of 
electromagnetic energy. 

In addition, the interaction between the jet against a disk wind 
or stellar envelope is also an unlikely source of current sheets 
(or numerous shells moving at varying relativistic speeds) due to 
toroidal-field-dominated relativistic jets tending to suppress bound- 
ary layer instabilities as discovered via 3D relativistic MHD sim- 
ulations of AGN jets (Keppens et al. 2009; Mignone et al. 2010), 
although specific studies for GRB jets should be performed. 



8 CONCLUSIONS 

We set out to explore a jet dissipation mechanism, denoted a "re- 
connection switch," that relies on the growing evidence that cur- 
rent sheets dissipate at different rates in collisional and collision- 
less plasmas. GRB jets naturally transition from being collisional 
to collisionless at large radii where dissipation is initiated near the 
jet photosphere due to a transition to fast collisionless reconnection. 

Our picture corresponds to a highly magnetized jet that is pre- 
sumed to be launched by a rotating neutron star or black hole. Cur- 
rent sheets are assumed to be generated by being injected into the 
jet by some dynamo action near the jet base or by instabilities at 
large radii. These reversals form a complex of numerous current 
sheets that dissipate slowly at small radii when the plasma is col- 
lisional and rapidly at large radii when the plasma becomes colli- 
sionless. The fast Petschek-like collisionless mode occurs once the 
ion skin depth is larger than the collisional current layer thickness, 
which allows the disruption of the thinner collisional layer geome- 
try and enables the fast Petschek-like reconnection geometry. 

For typical long-duration collapsar GRB parameters, recon- 



nection stays collisional until the jet transitions to fast collision- 
less reconnection at r ~ 10 13 -10 14 cm. Even if the reconnection 
switch mechanism fails to be valid and dissipation simply proceeds 
at the somewhat fast reconnection rate of v r < 0.01c, then, regard- 
less, all of our calculations remain intact because the dissipation 
radius (where significant electromagnetic dissipation occurs) is be- 
yond the expected transition radius for many models. However, if 
v r > 0.01c, then the reconnection switch is required to avoid com- 
plete photon thermalization in many models. 

Between the transition and dissipation radii, the Lorentz fac- 
tor is y ~ 100-1000 and the electromagnetic energy flux in the 
jet exceeds the kinetic energy flux by factors of 5-10. Assuming 
this energy is mostly dissipated and radiated away instead of ac- 
celerating the jet, then this allows the electromagnetic dissipation 
to produce a prompt GRB luminosity that is equal to (or even ex- 
ceeds) the associated afterglow emission, which could be what is 
observed but is difficult to obtain in the internal shock model (Will- 
ingale et al. 2007). Because y9j ~ 10-20, the afterglow can exhibit 
jet breaks. The jet is electromagnetically-dominated up to the tran- 
sition radius. However, a reverse shock can be present because dis- 
sipation decreases the electromagnetic energy flux to roughly the 
kinetic energy flux by the radius where an external shock forms. 

In our reconnection model, the number of prompt GRB pulses 
is suggested to be related to the number of current sheets embed- 
ded in the jet by the time the jet enters the dissipation regime and 
also related to the timescale for reconnection and radiative emis- 
sion to occur. For typical collapsar parameters, pulse durations are 
of order Is and range from O.OOls-lOs. Pulses have a rise time as- 
sociated with the transit time through a current layer as associated 
with heating, while the decay time is the timescale for dissipation 
of magnetic flux in the jet. These timescales are associated with 
a fixed transition/dissipation radius during the entire event so that 
pulses do not spread in time for fixed engine parameters. 

Significant electromagnetic dissipation only occurs once the 
flow has become optically thin for photons that reach the observer, 
and hence the compactness problem is naturally avoided. However, 
some portion of the dissipated energy could lead to some pho- 
tospheric emission for certain model parameters. This may help 
identify why some fraction of events seem to have thermal peaks 
(Ryde 2005). More work is required to obtain quasi-thermal emis- 
sion spectra (see., e.g., Giannios 2008). 

One goal of this work was to motivate future studies of recon- 
nection in the presence of (non-traditional) effects, including: rela- 
tivistic flow, pairs, radiation, super-critical electromagnetic fields, a 
complex of multiple current sheets, etc. This paper identifies impor- 
tant new directions in future reconnection research, and motivates 
them by identifying an new area of application of such studies to 
real astrophysical systems. 
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APPENDIX A: FULL JET STRUCTURE 

In this section, a generalized axisymmetric steady-state jet solution 
is presented as based upon approximate analytical solutions that 
have been confirmed numerically (Tchekhovskoy et al. 2008, 2009, 
2010b) using the HARM code (Gammie et al. 2003; Noble et al. 
2006; Mignone & McKinney 2007; Tchekhovskoy et al. 2007). 

In essence, highly magnetized jets are accelerated by the 
toroidal field gradient. The key to efficient acceleration is that field 
lines around the jet axis must redistribute so that an electromag- 
netic "exhaust nozzle" forms. At very small angles near the po- 
lar axis, the redistribution is highly non-self-similar, while slightly 
offset from the polar axis (where most of the energy flux is) the 
redistribution is approximately self-similar. In this approximately 
self-similar region, the toroidal field rapidly unwinds due to a drop 
in R 2 B p along field lines (Tchekhovskoy et al. 2009, 2010b), where 
B p is the poloidal field strength. One can obtain quite accurate yet 
simple analytical models of such jet solutions that apply to most of 
the jet undergoing acceleration containing most of the electromag- 
netic energy flux. The narrow non-self-similar region very close to 
the polar axis contains only a small fraction of the total energy flux, 
so only its effects on the rest of the mostly self-similar flow (not its 
own internal structure) need to be considered. 

Aspects of magnetized jets are well-modelled by a force-free 4 
(infinitely magnetized limit of ideal MHD) solution. The vector po- 
tential is roughly independent of rotation and given by 



RA 



<*,coll ' 



oc ( — I (1 - cos 9), 



where A^ fM is the vector potential for a collimating jet 5 , R = r sin 9 
is cylindrical radius, r fp is the foot point radius, < v < 1 is a 
free parameter, and 9 is the position angle in the domain. P coll is 
normalized by assuming that the radial field has the value B r (r fp ) at 
9 « n/2. The MHD invariants are constant on field lines (constant 
^coii), such that, e.g., the invariant field line rotation frequency can 
be expressed as Qf(r, 9) = Sl^(r( p , 9 lp ), where 9 fp is the angle of the 
field line at the foot point attached to the compact object or disk. 
For rapidly rotating black holes or neutron stars rf p f2p < 0.25c. 
Equation (Al) gives that the poloidal magnetic field obeys 



B r 

B r (r fp ) \r fp 



(A2) 



B r (r fp ) 



rtp 



tan((9/2), 



(A3) 



where the poloidal field strength is B p = ^B 2 r + Bj. Force balance 
considerations shows that the toroidal field is given by 



B,(r fp ) 



2r ip Clp 



tan(6»/2), 



(A4) 



which is consistent with the minimal torque condition that captures 
the effect of the Alfven critical surface (Narayan et al. 2007). 

The force-free Lorentz factor (y ff , such that y — > y ff for a 
pure force-free jet solution) follows the so-called first and second 



4 See, e.g., McKinney (2006a,c). 

5 This is accurate to 10% for force-free models with < v < 1.25 
(Tchekhovskoy et al. 2008). For ideal MHD models, it is accurate for the 
energetically dominant part of the jet (Tchekhovskoy et al. 2009, 2010b). 



asymptotic regimes and is defined by 
1 1 1 

v * V 2 + V 2 ' (A5) 

where in the first asymptotic regime 

7i - $ + (^yR/O 2 - (.ClFfpRff/c) 2 , (A6) 

where y is a free parameter to be determined later, and R fp = 
R(rf P , 9 fp ) with similar notation for other foot point quantities. The 
first asymptotic acceleration regime is dominated by the winding 
of the toroidal field. The second asymptotic regime has 
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where R c is the local poloidal radius of curvature of field lines and 
C ~ 3 (to order unity C depends upon the details of the solution). 
For large r, 
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(A8) 



The second asymptotic acceleration regime is dominated by the 
poloidal field curvature. The force-free Lorentz factor can be used 
together with MHD energy conservation constraints to obtain an 
accurate full MHD Lorentz factor. Given the total energy flux per 
unit mass flux (p, which is the upper limit of y) and the elec- 
tromagnetic energy flux per unit rest-mass flux (cr), then one can 
(Al) show that /i = y - ®tffl F B /(c 2 V47r) = y(l + cr) with the con- 
served quantity V47rO = B p /(p b u p ) = B r /(p b u r ) = B e /(p b u e ) = 
B^Hphiu^, - yRShfj). With the general estimate that y ff « fi/cr with 
cr — > as y — > /i, then one has that 

1 

Troll 

which accurately predicts the full Lorentz factor (i.e. y = y co ii) 
for collimating field geometries for MHD jets with v > 
(Tchekhovskoy et al. 2008, 2009). The value of y is chosen such 
that y = y fp at r = r fp . 

If the MHD jet was never confined or becomes unconfined, 
then the jet has v ~ and behaves qualitatively differently than a 
collimating MHD jet (Tchekhovskoy et al. 2010b). The first asymp- 
totic remains the same, while the monopolar solution has a modi- 
fied second asymptotic solution given by the cubic equation 

it/3 
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where Ci ~ 2, C2 ~ 0.4, 9 C ~ is the angle at which the jet passes 
the causality surface and y c ~ (p/ sin 2 9 C ) 1/3 is the Lorentz factor at 
the causality surface 6 (Tchekhovskoy et al. 2009). The maximum 
value of y 2 , m — » occurs as 9 — > 0, which leads to a maximally ef- 
ficient conversion of electromagnetic energy to kinetic energy. This 
behavior of the Lorentz factor near v ~ captures the crucial ef- 
fects of the fast critical surface and the related causality surface. A 
sharp transition 7 from the collimating to the monopole solution is 

6 As /j -> y fp ~ 1 these approximations for the causality surface introduce 
order unity errors. Generally [1 » 1 is assumed. 

7 As shown in Tchekhovskoy et al. (2010b); Komissarov et al. (2010), this 
sharp transition produces a rarefaction wave that generates the change in 
solution towards a monopolar one. A boundary layer at the outer angular 
edges of the jet undergoes a non-self-similar expansion, but this region con- 
tains only a small fraction of the total energy flux. 
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quite accurate treatment of the time-dependent ideal MHD solution 
(Tchekhovskoy et al. 2010b). The final Lorentz factor is then y co \\ 
inside r mono and otherwise follows the v = monopole solution 
with the second asymptotic given by y 2 , m . At very small 6 fp < 10~ 3 , 
these approximations can break down giving y < 1, in which case 
y = 1 is enforced because it is accurate for such regions. 

Now that the full Lorentz factor and magnetic field structure 
are defined, the ideal MHD invariants can be used to constrain the 
rest of the jet structure. The ^-component of the 4- velocity can be 
obtained from the conserved angular momentum flux per unit rest- 
mass fiux(/l = /Ifp = (Ru^-<bRB^)/(c y4n)) to obtain the conserved 
quantity if/ = if/^ = A - Clpfi = y - RflfU^/c 2 , where is the <f>- 
component of the 4- velocity thus given by 



y- ± 

RQ. ¥ 



(All) 



O is used to obtain u^ fp = Q. F ,f p R lp + B,,, Sp v pSp IB pSp . Foot point 
velocities are used to obtain yo(r fp , # fp ). 

The poloidal field-aligned component of the 4- velocity (u p = 



u 2 + u e) ls gi ven by the definition of the Lorentz factor using y 



2 _ 



such that u 2 



1 + (u 2 + up/c where v 



1/(1 - v 2 /c 2 ) 

c-V-D-4 

The rest-mass density can be obtained from the conservation 
of mass and magnetic flux such that 



Pb - Pi,tp 



A useful measure of the magnetization is the parameter 
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From the definition of fx one can show that p « yt p + 
4fc(n f>fp r fp /c) 2 sin(e f p)tan(e f p/2)/(y fp Up -fp ). For example, with 
f » Tfp ~ 1 one has that £ ~ 2p when 6 fp = n/2, Q, Ftfp = 0.25c/r fp , 
and 7f P ~ 1.15. In general, a foot point launching velocity of 
v P ,f P ~ c/2 is chosen (Tchekhovskoy et al. 2008), but it is unim- 
portant as long as 1 <, <k p. 

To enforce consistency with the ideal MHD invariants, B^ is 
recovered from the definition of p given the computed y and u^. The 
true solution's deviation of from Equation (A4) is typically less 
than 30%. Finally, the consistency with force-balance is checked 
by computing b 2 , which must be constant in 6 at each radius for 
a cold MHD jet to be in equilibrium. Generally, the region within 
7r/20 is found to be in approximately force balance. To obtain force 
balance across the entire jet, one iterates using the constancy of 
b 2 {ff) in the equation for p to obtain P co n (and so B r (6)) and y(6). 
Most of calculations in the paper focus on the most powerful part of 
the jet at large angles that is already in approximate force balance 
without this correction. This final solution satisfies all ideal MHD 
invariants exactly except the constancy of B^,/(p b {u^ - yRQ F )) once 
r > lOVp where this quantity has a relative error of less than 10%. 

The comoving electromagnetic pressure and energy density 
are p E M = "em = b 2 /(8n), where the comoving field is b p = (B p + 
(u ■ B/cJiu^/c^/y (ji index varies from t to <p). Using the fact that 
u ■ u = —c 2 and B' = gives the comoving electromagnetic field 
squared of b 2 = (B 2 + (u ■ B/c) 2 )/y 2 . For y » 1 and beyond the first 
asymptotic regime, this gives b 2 « B 2 /y 2 and b e b p « b r <<c b^,. 
So the field is dominated by the toroidal component in both the 
lab-frame and the comoving frame. 

The total power output of the jet is given by the integration 



over foot points, such that the single polar power output is 
C e 

Pj(0) = 2n SsmO'M'lLipb^Up], (A14) 
Jo 

with 6 allowed up to 6(9 fp = n/2), such that varying S rfp for 



£ » 1 leads to Pj cc 



r,fp 



for fixed other parameters. How- 



ever, because Pj and p are non-trivial (potentially non-monotonic) 
functions of free parameters, all plots will be shown with respect 
to the simple free parameters B rlp and At large distances be- 
yond the fist asymptotic, the electromagnetic jet power is p^ EM> = 

2n ^ r 2 dff sin 6' [b 2 y 2 v„] = 2n ^ r 2 dff sin ff [B 2 v p l 

We consider variations in f , S r fp , v, r mono , and 6 {p . The solution 
at large radii is insensitive to other parameters (e.g. fp.fp). 

These jet structure equations do not account for thermal en- 
ergy of GR effects. First, if the thermal energy is subdominant to 
the electromagnetic and rest-mass energy, then the above jet struc- 
ture equations are valid and remain accurate because only weak 
shocks occur in a highly magnetized flow. Second, GR effects ac- 
cumulate magnetic flux towards the black hole spin axis leading 
to a non-constant B rfp (6) (McKinney 2005b; Komissarov & McK- 
inney 2007; McKinney & Narayan 2007a,b; Tchekhovskoy et al. 
2010a). The spin enhancement of the magnetic field would also 
change how one would estimate B r j p from the mass accretion rate 
(Gammie et al. 2004; McKinney 2005b). However, these effects 
occur only within a few gravitational radii, and beyond this radius 
the flux deconcentrates without much difference to the non-rotating 
black hole case. This suggests a smaller true 9 fp for a rapidly ro- 
tating black hole with an accretion disk of height-to-radius ratio 
H/R -0.1 (as for a neutrino-dominated accretion disk) should give 
comparable results as generally choosing 6 fp ~ tt/2. 



APPENDIX B: DENSITIES, PRESSURES, AND EMISSION 
RATES 

In this section, the equation of state and emission rates for baryonic 
and radiative species are presented. For radiative species, the op- 
tically thick densities/pressures and optically thin energy/number 
loss rates are used with the two-stream radiative transfer slab ap- 
proximation presented in section B2. This approximation deter- 
mines the pressure and energy rates for general optical depths. 
These are then used in energy and force balance conditions that set 
the reconnection layer's physical geometry and reconnection rate 
in section C3. 

The distribution function for all species is assumed to be ther- 
mal except as modified by a radiative transport approximation to 
account for the range of optically thin and thick behaviors for radia- 
tive species. As tested using Equation (B8), the thermal assumption 
is often valid for the collisional layer. 



Bl Baryons 

A mixture of free nucleons (protons+neutrons) and a-particles can 
exist at the densities and temperatures of interest. Baryons are as- 
sumed to be at a single temperature T in thermal equilibrium. The 
baryons have a number density of n h = Pblm h , particle mass m b , 
and are assumed to be non-degenerate within the jet. Assuming nu- 
clear statistical equilibrium (NSE), the fraction of free nucleons is 



l,296p^Ce 1 ' 



(Bl) 
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(Woosley & Baron 1992; Kohri et al. 2005). The comoving 
timescale for establishing nuclear statistical equilibrium is dt^sE ~ 
p° b 2 exp(1.8 x 10" IT - 39)s in the electron degenerate regime and 
<^nse ~ [77(1.35 x 10 9 )]~ 5 s in the hot non-degenerate regime 
(Khokhlov 1989; Qian & Woosley 1996). These timescales are 
compared to the jet flow time in section 6.1. 

Given X nlic , the baryon internal energy density is 



"b = Pb 



3k^T I (3 - E m ( m )X nx + 1 



8m A 



1 



(B2) 



where £ rat i = £ , bin /(4ra iv c 2 ), the o--particle binding energy is 
E hm = 28.3MeV, m N c 2 « (m n + m„)c 2 /2 * 938.919MeV, and 
fe B is Boltzmann's constant. This directly accounts for the nuclear 
binding energy of cr-particles and the effect of photodisintegration 
directly within the equation of state. The baryonic pressure is then 



Pb = (T- IK, 



(B3) 



with T = 5/3, corresponding to a sum over all non-degenerate 
and non-relativistic baryons. The non-degeneracy assumption for 
baryons is ensuring by checking that n b <k (m b k B T/27tJi 2 ) 312 (Kohri 
& Mineshige 2002). Tests for the assumption that Y e = n p /n b ~ 1/2 
are discussed in section B5. 



B2 Two-Stream Radiation Approximation 

The energy densities and cooling rates for radiative species are ob- 
tained via the two-stream approximation limit of the Boltzmann 
equation (Rybicki & Lightman 1979; Hubeny 1990; Popham & 
Narayan 1995; Sawyer 2003). The theory of radiative reconnect- 
ing current layers is quite analogous to the radiative transfer prob- 
lem in accretion disk theory (Uzdensky & McKinney 2011), for 
which this approximation has been used extensively for both pho- 
tons (Popham & Narayan 1995) and neutrinos (Popham et al. 1999; 
Di Matteo et al. 2002; Kawanaka & Mineshige 2007). In this one- 
zone approximation, radiation is either lost from the slab or remains 
and is not redistributed within the slab. 

The two-stream approximation requires optical depths for ab- 
sorption (T abs ), scattering (r sca ), and their total 



Ttol 



(B4) 



For scattering species A interacting with a bath of scatterers 
species B with number density n B , scatterer density scale-height 
H B , and scatter cross section cr A _ BiSca , the scattering optical depth is 



T A,sca - / i n A-B,eflHB cr A-B,scn- 



(B5) 



Order unity changes in the momentum of species A are re- 
quired to achieve an effective collision corresponding to large-angle 
scattering and an effective exchange of energy. The actual number 
density of scatterers B given by rig is effectively reduced by diffu- 
sive (rather than direct) scattering, such that 



and 



tiA-Btfs = n B X[p A ], 



■ | Ap A Ap A 

X[p A ] ~ mm , 

I Pa \ Pa 



(B6) 



(B7) 



where p A is the momentum of species A. The center-of-momentum 
frame change in momenta is approximated as Ap A ~ mm(p A ,p B ), 
where p B are the momentum of species B. The min conditional 
in accounts for the fact that effective collisions occur when 

(Ap A )/p A > 1, while the momentum undergoes a random walk for 



(Ap A )/p A < 1. This correction is not applied to Coulomb scattering 
as computed later because this effect is already included directly in 
the cross section. 

Thermal equilibrium between a pair of particles is achieved if 
the transit time of plasma through the current layer is longer than 
their collisional energy-exchange time. This condition is written as 

L ° " 1 (B8) 



Cthermalization,A-B ~ v ec,A-B 1, 
Va 



where the effective collision energy-exchange rate due to order 
unity changes in the momentum of species A is 

VecA-B ~ V acA - B X[p A ], (B9) 

where is given by Equation (B7) and v acA - B is the actual col- 

lisional frequency for particle type A to collide with particle type 
B. In the non-relativistic limit, v ecf - e ~ n e T'g' :2 m l } 12 ~ v e c,e-i ~ 



{milm e ) 



1/2,. 



(m i /m ( ,)v ( , ci _ ( ,, where ions are scattered diffu- 



sively by electrons. The value of C t hermaiization,A-B is checked for all 
models in order to ensure sufficient thermalization, although even 
in the collisionless regime one expects a dominant thermal compo- 
nent (Giannios & Spitkovsky 2009). 

The absorption opacities are obtained using Kirchhoff's law, 
such that for each absorber with number density n B and some ab- 
sorption scale-height H B one obtains an energy and number density 
loss rates based absorption optical depth of 

Qb-a,oH b 



r u.A,abs = ^ n B H B CT A_s,„,abs - ^ 



1"n,A,abs = / . n B H B (T A-B,„,abs : 

V V W ftA«0,A 



z 



VhA u 0,A 

Rb-a,qHb 



(B10) 



(Bll) 



respectively, where u ,n are the optically thick (i.e., those corre- 
sponding to the thermodynamic equilibrium at a given temperature) 
energy and number densities, respectively, and v fs is the free-stream 
velocity (e.g. v fs = c for photons and neutrinos). (Notice that the 
subscript of only indicates the quantity is some known quantity 
used to construct new quantities that are valid at general optical 
depths.) The cross section cr A _ B :ibs is for the inverse reaction to that 
operating at the optically thin rate Qb-a,o, Rb-a>)- 

Now that the absorption and scattering optical depths have 
been obtained, the two-stream approximation can be used to ob- 
tain densities and rates valid for general optical depths. The internal 
energy and number density are given by 

u = u g[r u ], (B12) 
n = n g[T„]. (B13) 

with 



t«/2+1/V3 



(B14) 



T tot /2+l/V3+l/(3r abs )' 
and energy and number surface fluxes lost from the slab of 

F u = u fs M /![r„], (B15) 
F„ = v fs n h[T n ], (B16) 

with 

1/3 



h[T] 



(B17) 



T tot /2+l/V3 + l/(3r abs ) 
The volumetric energy and number density rates are then given by 



F„ 



«,abs 



(B18) 
(B19) 
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respectively, where the absorption depth is 

T abs 



^absl^abs] = 



(B20) 



dT dbs /ds 

where from Kirchhoff's law one obtains dT u ^ s /ds = 
Zs «flO"B,«,abs = 2s 2s,o/(ffc"o) = 6o/(«ft"o) and 
dT„# bs /ds = 2s«B cr B,n,abs = 2s ^B,o/(Vf s «o) = Ro/(v ts n ). 
Notice that g, W become the optically thin energy and density rates 
go^o when t <k 1. 

B3 Optical Depth in a Relativistic Jet 

In order to compute the optical depth in the prior section, one re- 
quires some estimate of the density scale height for scatters, given 
by H, that is applicable for a relativistic jet. 
In the lab-frame the optical depth is 

t = J cdt'[crn] = J dr [crny( 1 - p cos 0)], (B21) 

where cdf = dr' is the comoving time difference, dr' = y(l - 
P cos B)dr is the lab-frame radial difference, n is the comoving num- 
ber density of particles, y is the Lorentz factor, cf! is the 3-velocity, 
and 8 is the angle in the lab-frame. An angle-averaged invariant 
cross section (<x) is assumed to be used. For a fluid with y » 1 in 
the lab frame, a typical photon emitted isotropically in the comov- 
ing frame is emitted parallel to the jet in the lab-frame giving 



f 



dr[ncr/(2y)], 



(B22) 



where Ar is the lab-frame distance the photon traverses through 
the medium of density n. For photons to escape a radiating slab of 
comoving size L , Ar ~ 2yL per layer would be chosen. 

In our case, the dissipating complex of current sheets consists 
of numerous narrow dissipation regions covered in a photon-pair- 
neutrino photosphere (when a species is optically thick). As shown 
in section 5.4, the range over which the electromagnetic energy is 
dissipated into these species is given by Ar ~ cyA /(2v r ) for all 
m modes and some I modes. Assume that radiation densities drop 
as roughly r~ 2 (accurate for photons up to relevant radii), and as- 
sume that the Lorentz factor and cr(p, T) vary as some power of r 
(including constant with radius). The opacity integral then gives 



Trad 

crn md 



ysin(0,) 2v r , , 



Ri 



cA 



(B23) 



where n mi is the number density of radiation at radius r. For down- 
stream baryonic-associated electrons, a similar calculation gives 



/ Tbaryonic \ ' ysin(flj) 
\ Cr^baryonic / Rj 



(R'jY 



(B24) 



where baryons are assumed to extend to large radii such that Ar — > 

CO. 

So we have determined that H ~ A' for scatterers of radiative 
species generating by the current layers, and H ~ R'. for scatterers 
consisting of baryonic-associated electrons downstream in the jet. 
A more detailed relativistic radiative transfer calculation is left for 
future work (see, e.g., Meszaros 1992; Beloborodov 2010b). 

The two-stream approximation can now be applied to photons, 
pairs, and neutrinos in order to determine their behavior in the op- 
tically thick and marginally optically thin regime when the single 
temperature approximation is still good to order unity. This two- 
stream approximation determines the densities (w), pressures (p), 
energy rates (Q), and number rates (R) for photons, pairs, and neu- 
trinos as obtained from the base optically thin rates (Qo,Ro), base 



optically thick densities (w , «o), and scattering cross sections (cr sca ) 
written down for each species in the next section. 

B4 Photons 

Photons are created within the dissipative current layer by pro- 
cesses such as bremsstrahlung (e.g. free-free emission, free-bound 
emission), cyclo-synchrotron, pair annihilation, radiative pair an- 
nihilation, and double Compton scattering. Most photons are cre- 
ated inside the dissipative current layers, and photons can travel 
across field lines and interact with other layers. For photons, the 
two-stream approximation requires knowing the scattering and ab- 
sorption opacities as computed below, and one also needs the opti- 
cally thick limit for the energy and number densities. The optically 
thick photon internal energy density is 



«0o 



: a y T\ 



(B25) 



where a y = n 2 k^/(l5(tic) 3 ). The photon pressure is then given by 
Po, y = «o, r /3 as for a T = 4/3 ideal gas for adiabatic constant T. 
The optically thick number density of photons is 

n , y = (2<r[3]/7r 2 )(V£ B )r 3 ~ (1/4)0W(W). (B26) 



B4.1 Photon Scattering 

The scattering opacity for photons is 



<T e s(« r -pairs,effAo + n y _ efiS R'j), 



(B27) 



for an electron scattering opacity <r es , where n r _ pa i rs ,i;fr is the ef- 
fective number density of pairs and n y - efi g is the effective num- 
ber density of electrons as computed from Equation (B6), and the 
scale-heights for scatterers (A' and /?'.) are given by Equation (B23) 
and Equation (B24). Radiative effects discussed later can lead to an 
enhancement of the baryon and associated electron density, but the 
opacity is an integral of density that is roughly fixed for a conserved 
amount of mass across the jet. 

In the sub-critical QED regime, the thermal spectrally- 
averaged cross section for the Klein-Nishina effect is roughly 



o~kn \ 



1 + 



3m e c 
Sk B T 



(B28) 



where cr T is the Thomson cross section. 

In the super-critical QED regime, the extraordinary mode (E- 
mode, electric polarization perpendicular to the magnetic field), 
electron scattering cross section depends upon whether the particle 
energy is above/below the rest-mass of electrons and above/below 
the energy of the first Landau level given by hw Be (l) = ((m e c 2 ) 2 + 
2hcq\b\) 112 - m e c 2 (Silantev & Iakovlev 1980; Meszaros 1992; 
Thompson & Duncan 1995; Lai 2001). For the E-mode, the thermal 
spectrally-averaged cross section is 



^| -1 + 



^f5imi e c 2 k B T 
hcqb 



(B29) 



which only applies for super-critical field strengths. The high-field 
suppression of the scattering cross section for the E-mode and the 
efficient conversion of the ordinary mode (O-mode) to E-mode 
(Meszaros 1992; Thompson & Duncan 1995) means that radia- 
tive emission is dominated by the E-mode with this scattering cross 
section when above the critical field strength. Photon splitting and 
merging are assumed to be in detailed balance in the optically thick 
regime (Thompson & Duncan 1995), and often the thermal photons 
are optically thick for super-critical fields. Defining an interpolation 
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parameter e cB = exp ((-m e c 2 )/(ha> Be (l))), the total electron scatter- 
ing cross section for any T and b is 



cr B e eB + cr KN (l - e eB ). 



B4.2 Photon Absorption and Emission 



(B30) 



The absorption opacity is determined from the optically thin emis- 
sion rate and Kirchhoff's law. Because the true number density of 
photons is not required to be accurate in this work, r y abs « r„ r abs ss 
T«, r ,abs is set an d only energy density loss rates are considered. Ther- 
mal free-free, thermal synchrotron, and thermal pair annihilation 
are considered. The total photon absorption optical depth is 

T r ,abs ~ (fff + 0~ synch)(Wpair S Ao + n e R'j) + CT e + e - ^ yr n p . dirs A' Q , (B31) 

as due to, respectively, free-free, synchrotron, and pair annihila- 
tion given below. While the non-pair-producing electrons (with 
number density n e ) outside the current layer's absorption photo- 
sphere may be initially cold in the jet, those electrons are as- 
sumed to be heated by a sufficient number of photons bringing 
that portion of the jet into thermal equilibrium if the region is 
optically thick. The local differential absorption optical depth is 

dTy^bslds = (0" ff + <T sync h)(Mpairs + «e) + CT e+f-^yyWpairs- 



B4.3 Free-Free 

The optically thin free-free emission rate that includes relativis- 
tic effects, electron/positron-ion collisions, and electron/positron- 
electron/positron collisions is 

\2 



<2i, ff ~ 1.4x10-^7- 



-27 T 1/2 



1 + 



KE e 



(B32) 



where KE f « 3k B T is the average kinetic energy of relativistic elec- 
trons/pairs, n e ~ n p ~ n b is assumed, and protons are assumed to 
be non-relativistic (Bottcher et al. 1999). Kirchhoff's law gives 

Giff 

^o, ff = — , (B33) 

«e,totCM0, r 

for the cross section per emitting-absorbing particle. In the super- 
critical QED regime, the thermal spectrally-averaged cross section 
is suppressed such that 

_ N -2 

Off 
\°~0,B ) 



I o-ff Y 1 ^f5nm e c 2 k B T \ * 

Wofi) { ncqb } 



(B34) 



(Thompson & Duncan 1995; Lai 2001). Note that one can define 
the QED optically thin emission rate, Q off , from Kirchhoff's law: 
o~g = 2o,ff/(«e,totC«o,y)- These equations give, as required, a free- 
free contribution to Q g of g ,ff in the optically thin limit. 

B4.4 Cyclo-Synchrotron 

The synchrotron (and approximate cyclotron) energy density loss 
rate integrated over all angles, frequencies, and over an isotropic 
distribution of thermal particles with Lorentz factor y e and pitch 
angle 6 is 

J r>7i pes 
sine | dy e [f(y e )j(y e ,6)], (B35) 
9=0 Jy e =l 

where the frequency-integrated emissivity (erg/s) for each pitch an- 
gle is 

2q 4 b 2 (y 2 - l)sin 2 <9 



3m 2 c 3 



(B36) 



and for simplicity the non-degenerate relativistic electron thermal 
distribution, 

1 



f(7e) ■ 



ylfieexpi-yJ&X 



(B37) 



0K 2 (l/0) 

is used, where ® = k B T/(m e c 2 ), B e = -^1 - 1/y 2 , and K is the 
BesselK function (Rybicki & Lightman 1979; Mahadevan et al. 
1996; Ozel et al. 2000). In cases when the gas is degenerate, syn- 
chrotron is found not to be crucial to this study. The electron and 
photon azimuthal angles have both already been integrated over. 
Kirchhoff's law gives 



0~ O.synch 



Q 



1 ,synch 



(B38) 



n em cua, y 

QED effects are important for |&|/Z?qed ^ 0.01 where &qed = 
{m e c 2 ) 2 1 {cqh) and always lead to a suppression of cr 0synch . For 
y e — 1 > |b|/fcQ ED , well-defined expressions exist for computing a 
QED version of <x synC h (e.g. consider the complete integral of equa- 
tion (31) in Baring 1988). However, the full integrals are compu- 
tationally expensive. So instead of directly using the full QED ex- 
pression, we numerically derived a suppression factor fitting func- 
tion. The suppressed cross section is 

I (m e c ) I 
1 + 1 _ 1 , (B39) 



'synch 



(m e c 2 Y 



\ Co.synch / V ^fiixhcqbk B T ) 

which is found to be accurate to order unity for |£>|/Z>qed ^ 0.01 and 
> 0.01 as is sufficient for this study. The asymptotic suppression 
factor is cc (r|Z?|)~ 4/3 , such that both high field strengths and high 
temperatures induce a suppression effect. Note that one can define 
the QED optically thin emission rate, Qxsyndu from Kirchhoff's 
law: tr synch = 2o,synch/("e,totC«o, r ). These equations give a contri- 
bution to Q g of Qo.synch in the optically thin limit. 

The total inverse Compton power is roughly a factor 
(« y /w EM )(<T es /cr T ) of the synchrotron power. In this work, u y < « EM 
and cr es < o~ x, so that the total inverse Compton power is typically 
weaker than the total synchrotron power. In the transition from 
optically thick to optically thin radiation, synchrotron at low fre- 
quencies can be self-absorbed leading to a dominant Comptonized 
emission (Giannios 2008; Lazzati & Begelman 2010). This effect, 
a calculation of (Comptonized) spectra, and related QED effects 
(Harding & Lai 2006) are left for future work. 

B4.5 Pair Annihilation into Photons 

Electron-positron pair annihilation into photons for thermal pairs 
has a number density rate of 



Ro,e+e-^ry ~ o^TCH+n- 1 + 



20 2 



(B40) 



ln(1.120+ 1.3)) 

for each thermal pair unit (1 electron and 1 positron), where = 

k B T /(m e c 2 ),n + = (n em -n e )/2 = n pairs /2, and = (n efot +n e )/2 = 
Wpairs/2 + n,, (Svensson 1982; Sikora & Zbyszewska 1986). The cor- 
responding energy density loss rate is 

Qo,e+e-^yy ~ £p a i IS Z?o,e+ e~~ -+77' (B41) 

where ,Ep a j rs ss 2(m e c 2 + u e /n e ), and u e is the electron/positron in- 
ternal energy density given later. By Kirchhoff's law, 

QED effects on the optically thin pair annihilation rates are ne- 
glected, but only tend to be important at very low baryon loading 
of the jet. 
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B5 Electrons and Positrons 



B 5. 1 Electrons and Positrons EOS with QED Corrections 



Electron-positron pairs can be created by processes such as pho- 
ton annihilation, one-photon pair production in super-critical fields, 
ion-electron collisions, electron-electron collisions, ion-photon 
collisions, and electron-photon collisions. The pairs do not read- 
ily cross field lines and are stuck within the dissipation region until 
they can annihilate into photons and spread throughout the photon 
photosphere or travel along field lines within the layer. 

In thermal equilibrium, the asymmetry between electrons and 
positrons is represented by the electron chemical potential p e = 
p e - = —n e *, which is determined by the condition of charge 
neutrality among protons, electrons and positrons given by n e = 
n e - — n e + = n p , where n p is the total number density of protons 
(free or bound). The degeneracy parameter of electrons is given by 
t] e = jiJk^T, where tj e > 1 implies electrons are degenerate. A 
value of Y e = n e /n b = 1/2 is assumed, but the /^-equilibrium value 
of Y e is computed from l/(l/Y e - 1) = n p /n„ = exp(2 - n e ) with 
Q = (m n - m p )c 2 /(k B T), which assumes only free nucleons are 
present (Kohri & Mineshige 2002). This gives a test of whether the 
assumption of Y e = 1/2 is violated. 

The Fermi-Dirac distribution function of particles in thermal 
equilibrium is 



ME) 



1 



(B43) 



e (£/(* B ?>i;.t) + 1 ' 
The number densities of electrons and positrons are 

i r°° 

"< ± = nHhc7 I MP (B44) 

where p = pc and E e = ^jp 2 + (m e c 2 ) 2 . For the total number den- 
sity of (free and bound) protons (n p = Y e p b /m b ) and T, the value 
of p e is iteratively solved for from the definition: n e (p e )jn b = Y e . 
Then n e - and n e + and any other quantities where p e appears can be 
computed. The electron-positron pressure is 



Pe* 



1 



3n 2 (hc) : 



ife<E e )) 



(B45) 



The electron rest-mass plus internal energy gives an internal energy 
density of 

i r°° 

"* ± = ^7iT^ dp[p 2 E e (fe±(E e ))]. (B46) 
n 2 (hcf Jo L J 

The non-degenerate limit gives 



p e = n e k B T, 



and 



r- 1 



(B47) 



(B48) 



where the parameter = (k B T) / (m e c 2 ) is used to linearly interpo- 
late from r = 5/3tor = 4/3 such that at and beyond a value of 
unity this gives T = 4/3. 

Then, the pair pressure, internal energy density, and number 
density are 



PO.pairs = Pe- + Pe 
M 0,pairs = u e- + M e+ 
«0,pairs = n e* + «e" 



~Pe 
■ U„, 



2n e 



(B49) 
(B50) 
(B51) 



For more details, see Kohri et al. (2005). Note that the sum of, e.g., 
Pe + Ppmrs that enters the total gas pressure is unaffected by this 
decomposition. 



In the super-critical field regime, the pair internal energy density 
depends upon whether the particle energy is above/below the rest- 
mass of electrons and above/below the energy of the first Landau 
level. Let ei = exp ((—m e c 2 )l(k B T)) and introduce an interpolation 
factor e 2 = exp ((-k B T)/(haj Be (l))) (Meszaros 1992; Thompson & 
Duncan 1995). Then, the following approximate QED expressions 
are used if e 2 > 0.1, and otherwise the more accurate non-QED 
expressions given previously are used. 

When hoj B e(l) •« k B T and k%T » m e c 2 the pairs behave like 
classical radiation with an internal energy density 



"OApairs ~ (7/4)Wo, 7 ei(l - e 2 ), 



(B52) 



which gives a number density of pairs roughly equal to the number 
density of m e c 2 photons. When hco Be (l) <<c k E T and k B T <K m e c 2 , 
then the pairs are classical and non-relativistic with an internal en- 
ergy density 



2 l < 2 (m e c 2 f (k B T\ 312 

"O.g.pairs ~ | 9 | Gl(>- ~ Vl)- 



(Hey 



(B53) 



When hcL) Be {i) » knT and k%T » m e c 2 , then photons and quan- 
tized relativistic pairs are in detailed balance with an internal en- 
ergy density 

1 ficqb 



"O.Cpairs 



12 (he) 



■(k B T) eie 2 . 



(B54) 



When nu>Be(l) » k^T and k^T <K m e c 2 , then the pairs are quan- 
tized and non-relativistic with an internal energy density 



w 0,D,pairs 



(hcqb)(m e c 2 ) 2 



1/2 



eie 2 . 



(B55) 



(2^ 3 ) 1 / 2 (fc) 3 

Across these four regimes, a sufficient interpolation procedure is to 
simply sum all terms together to obtain an internal energy density 



"O.pairs ~ "0,A,pairs + "O.B.pairs + "O.Cpairs + WO.D.pairsi 



(B56) 



pressure 

Po.pairs ~ "0/l,pairs/3 + 2w ,S,pairs /3 + "0,C,pairs/3 + 2w ,O,pairs /3, (B57) 

and pair number density 

/'O.pairs 



'^0,pairs 



k B T 



(B58) 



These expressions are only accurate to order unity for k B T > m e c 2 
and otherwise less accurate, but the QED effects tend to only be 
relevant at high temperatures. 

B5.2 Opacity Effects on Electrons- Positrons 

The pair production processes yy — > e + <T, ye — > ee + e~~ , and ee — > 
eee + e~ are considered, where the ep — > epe + e~ and yp — » pe + e~ 
have been shown to be less efficient for thermal plasmas considered 
in this paper (Zdziarski 1982; Sikora & Zbyszewska 1986). 

Pairs act as radiation when annihilation into photons domi- 
nates pair creation and pairs fill-in the region above the current 
layer with a density scale-height of order Eq as for photons. In the 
limit that Kirchhoff's law applies, the creation of photons occurs 
via the inverse reactions to yy — > e + e~ and to ye — » ee + e~. As 
seen below, these dominate the inverse of ee — > eee^ e~ except at 
k b T » m e c 2 , such that any absorption opacity leads to photons that 
can readily cross field lines and come into equilibrium throughout 
the photon photosphere of size L . Pairs also act as radiation (with 
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the addition of Coulomb interactions) when pairs flow down the 
current layer along straight field lines that open-up after a length 
L . Then, both across and along the layer, the photon and pair den- 
sity scale-heights are order L . In the collisionless limit, the pairs 
are determined by their prior history instead of equilibrium, but this 
work only seeks to find the collisional layer structure. 

The two-stream approximation is applied to pairs whether they 
annihilate and travel above the layer or they travel down field lines 
along the layer. Order unity factors in equations (B14,B17) would 
slightly change depending upon the allowed trajectories, but such 
minor changes are ignored. Kirchhoff's law is used as usual with 
angle and energy averaged rates given below, where recall that 
{crv)/c = (1 + Sn)R n /(cnin 2 ) for species 1 and 2 and number 
density rate R [2 (see., e.g., Weaver 1976). For simplicity, the un- 
derlying photon, electron, and pair distributions are assumed to be 
thermal, although the photon densities are modified by the two- 
stream approximation to account for optical depth effects. Because 
the (fully or marginally) optically thin energy density and pressure 
of pairs are not required to be accurate in this work, while the num- 
ber density of pairs in the marginally optically thin limit should be 
somewhat accurate, the number rates are treated directly while the 
energy rates are approximated. 

Consider the yy — > e + e~ process. Following the notation 
and assumptions of Weaver (1976), let C = exp(-fi y /(k B T)) with 
photon chemical potential p. y , <t> = m e c 2 /E y for the center-of- 
momentum photon energy E y , = k B T /(m e c 2 ), and x = 1/(0®) 
such that xk B T = E y , then 



3o"r 



<D 2 (2 + 2<D 2 - O 4 ) coslT 



<f>' 



(B59) 



for <D < 1 and cr yy ^ e - -> for O > 1, where <D' = (1 + $ 2 )(1 - 
$ 2 ) 1/2 and cosh - ' is the inverse and not reciprocal. Then the num- 
ber density rate for each pair unit is 



R ,ry^e- = n 2 y c ]T (V^/C" +1 )-'x 



(B60) 



f dx\ 

Jx=0 



dx x 4 o~ v 



where K is the BesselK function (Weaver 1976), and fi y — > such 
that C = 1 is enforced because the variation in photon number 
density is subsumed into the prefactor number density of photons 
n y (which would otherwise have been n 0y and one would have to 
solve for ji y instead of using the two-stream approximation). The 
energy density loss rate is 



Go,, 



E VV R, 



yyR-Q,yy—>e + e 



(B61) 



where the unit of pairs has thermal energy of roughly where E y 
2E y and E y ~ u y /n y . 



The ye —> ee + e process has a fitted cross section 



- A = l(T 3 (>'-4) 2 x 



[5.6 + 20.4(y - 4) - 10.9(y - 4) 2 
- 3.6(y - 4) 3 + 7. My - 4) 4 ], 
r ye ^ee+e-,B = 0.582814 - 0.29842? + 0.04354/ 
-0.0012977/, 

3.1247 - 1.3397? + 0.14612/ 



Vye-*ee + e ,C 



Cye^>ee + e ,D 



1 + 0.4648y + 0.016683/ ' 
; (841n(2y)-218)/27 
-1.333 \n 3 (2y) + 3.863 ln 2 (2y) - 11 ln(2y) + 27.9 



(B62) 

(B63) 

(B64) 
(B65) 



where y = y r = E y /(m e c 2 ) is the variable photon energy per elec- 
tron rest-mass energy, such that 



3crj-a 



8/r 



if y < 4 
-.a if 4 < y < 4.6 



<Tye^ee+e-,B if 4.6 < >' < 6 
tT' re -,«+ e -,C if 6 < >' < 14 

cr ye ^, ee + e -j} otherwise 



(B66) 



where a is the fine structure constant (Stepney & Guilbert 1983). 
The underlying photon distribution is assumed to be Bose-Einstein, 
such that the number density of photons is 



2(k B T) 



{he. 



' n=l 



with distribution 



fy 



(Pycf 



(ftc) 3 «o, r Ce<PrO/<t B r) _ j ' 



(B67) 



(B68) 



for variable photon momentum p y = E y /c. The number density 
rate for each pair unit is then 



^0,ye—>ee + e ^e. 



2K 2 [l/0] 



(B69) 



f dydy r I 

Jv=0,v r =4 L 



dydyAy f y y r cr ye ^ e , 



y/yy+yr/y "' 

[y r ]e ^~ 



Jy=0,y r =4 

(Sikora & Zbyszewska 1986; Zdziarski 1982; Svensson 1984, 
1987). The energy density loss rate is 



Qo.y 



e - % E ye Ro_ 



(B70) 



where E ye ~ E y + KE e and KE C ~ u e /n e . 

The ee — > eee + e~ process has a number density rate for each 
pair unit of 

= 8.4 x lO" 6 n 2 lot cr r c(ln(0)) 3 , (B71) 



R a 



for © > 1 and R 0tee ^ eee + e - — > for < 1 (Sikora & Zbyszewska 
1986; Svensson 1984, 1987). The energy density loss rate is 



Go. 



EpaRt 



ee 1 ^0,ee-*eee + e 



(B72) 



where E ee ss 2KE C . 

The total optically thin energy and number density rates are 



^O.pairs — ^0 



Go.pairs = Qo,yy^e+e- + Qo.ye^ee+e- + Qo.ee^eee+e 



(B73) 

(B74) 
(B75) 



The scattering optical depth is 

^"pairs.sca ~ ^"es^pairs-y.efF^O ^"c(^pairs^o ^eRj\ 

where n P air S -y,eff is the effective number density of photons com- 
puted from Equation (B6) using p y ~ (l/c)(u y /n y ), and <t c is the 
Coulomb scattering cross section given later by Equation (C2). The 
absorption optical depths are given by Kirchhoff's law as 

2^0,pairsAQ 



' rt.pairs.abs 



^"w.pairs.abs ~ 



WfsWO.pairs 
Go.pairsAg 
^fs ^O.pairs 



(B76) 
(B77) 



where the factor of 2R 0pilks appears because /?o,pairs corresponds to 
the rate to produce a pair unit, while « ,pairs is the total number of 
pairs separately counting electrons and positrons. H a \, s = AJ, as hap- 
pens because there is only one depth involved for pairs and pho- 
tons involved in the absorption depths. The effective free-stream 
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velocity (w fs ) for pairs traversing across field lines over the scale- 
height A' Q is determined by the distance per unit time traveled as 
photons and pairs. The fractional distance traveled as photons is 
L y ~ c/(R 0i77 ^ e + e - + Ro, y e^ee+e ) and the fractional distance trav- 



eled as pairs is L pii i 



The total fractional travel 



time is T ~ l/R 0y e+ e -->yy + l/(fio i7rw t r + Ro, re ^ ee + e -), such that 
the effective average free-stream velocity is t> fs ~ (L y + L pairs )/7\ 
Across the field lines v e ~ due to particle gyrations around mag- 
netic field lines, while parallel to the field lines, v e ~ c -y/1 - Tfy? 
is the electron-positron thermal speed. In addition, along the length 
of the layer, the pairs can be adverted at the speed v A ~ c. For 
the temperatures and densities considered, a good approximation is 
found to be vt s ~ v e ~ c, where lower v e correspond to low temper- 
atures where pairs are not dynamically important. The total opacity 
is Tpairs.tot = Tpairs,sca + Vi rs ,abs- The two-stream approximation gives 



Ppairs — /'(Xpairsgl^w.pairsL 
^pairs — ^O.pairsglTw.pairsL 
^pairs — ^0,pairsglTtt,pairsL 
Gpairs — ^fs w O,pairs^[Tw,pairsL 



Rn 



«fs«0,p 



s /z[T tt n a i rs ]. 



(B78) 
(B79) 
(B80) 
(B81) 
(B82) 



QED effects (e.g. 1-photon and modifications to 2-photon an- 
nihilation) are only important when the pairs are optically thick to 
both absorption and scattering, so the QED effects (see., e.g., Bar- 
ing 1988; Baring & Harding 1992) on the optically thin rates need 
not be considered for the density of pairs. 

B6 Neutrinos 

At the highest densities (p b <; 10 10 g/cc) and temperatures (T > 
10 10 K) considered, neutrino emission is the dominant source of 
cooling and pressure. Kirchhoff 's law is used as usual, except be- 
cause the true number density of neutrinos is not required to be 
accurate in this work, T v>abs ~ t„ y abs ss r M>v>abs is set and only energy 
density loss rates are considered. 

We follow Kohri & Mineshige (2002) (see also Kohri et al. 
2005), except their energy density loss rates assume that there ex- 
ists an optically thick thermalized photon and pair bath, while in 
this work the photon and pair densities are reduced by a factor g[T r ] 
and g[T pails ], respectively. When photons and pairs are involved in 
the reaction, a g[r] factor is applied for that given number density 
(oc T 3 ) as it enters the original integral, so that these rates are con- 
sistent with the number density of photons and pairs. Without this 
correction, then (for example) the pair annihilation rate would be 
erroneously large at high temperatures in the regime where pho- 
tons and pairs are optically thin. This assumes, as accurate in this 
work, that pairs dominate electrons in number density when neutri- 
nos are being produced. The neutrino energy density loss rate for 
capture of non-degenerate pairs on nucleons is then 

< 9.2 x W 33 T 6 n (fi w X mc )g[r p ^ s l (B83) 



and on degenerate pairs is 

Qw„ * l.i x lo 3 ^^^], 



(B84) 



where Ne — > denotes a sum of processes e + + n — > p + v e and 
e~ + p — > n + v e and — > Ne is used to denote the sum of their inverse 
reactions. The value of rj e is used to linearly interpolate between 
regimes, such that at rj e = 1 only the degeneracy expression is used. 
Neutrino pair production by annihilating pairs gives 

n 2 



Go,e 



+e —>v+v 



■ 4.8 X lO^rj'.gtTpairs] 2 



and is negligible in the electron-degeneracy regime. Neu- 
trino pair production by non-degenerate free nucleon-nucleon 
bremsstrahlung gives 

Go 1.5 x 10 }i T™(p l3 X mc ) 2 . (B86) 

The plasmon process gives 

Go,^v + v - 1.5 x I0 32 r» r *e-»(1 + y P )^2+ y|^-jg[Tp (df ,]g[T r ], 

(B87) 

where y is a photon interacting with electrons and y p = 5.565 x 
10~ 2 [(7r 2 + 3;7 2 )/3] 1/2 . Then the total optically thin neutrino energy 
density loss rate is 



Qo.v - Qo,Ne^ + Qo,, 



e + +e — >v+v 



+ Qojl+n->, 



Q,n+n— >n+n-hv+v 



+ Qo.y- 



(B88) 



The neutrinos can be optically thick at sufficiently high den- 
sities/temperatures, which is important to include because an arti- 
ficially high neutrino cooling rate would spuriously lead to com- 
pressible solutions for the collisional layer (see section C3). The 
electron type neutrino is treated most accurately among all the neu- 
trino species because it generally dominates the energy density loss 
rates. Both anti-neutrino and neutrinos are treated using a single 
opacity and all neutrino chemical potentials are assumed to be zero. 
The scattering optical depth is 



Ty.sca - °" v,bt( n v-bt,efsWj + V ve* (Mv-e±,i;flA + n e R'j), 



(B89) 



where n v _bf,eff is the effective number density of free baryons, 
«bf = X mc n b is the number density of free baryons, and n v -e±,eff is 
the effective number density of electrons as computed from Equa- 
tion (B6). Free nucleon scattering has 



o- Vjbf = 7.7 x lO-"m b (C s „Y p + C s , n Y n )T\ 



(B90) 



where Y„ = Y e and Y n = 1 - Y„, C StP = [4(C V - l) 2 + 5a 2 ]/24 and 
C V1 = (1 + 5a 2 )/24, with vector coupling C v = 1/2 + 2s'm 2 6 w , 
a a ~ 1.25, and Weinberg angle is sin 2 6V = 0.23. For electron- 
positron scattering 

(B91) 



<T ve ± = cr (l + 



1 



(C v + C A f + -(C v - C A ? 



k B T , 

m.c 1 



where C A = 1/2 for electron neutrinos and C A = -1/2 for elec- 
tron anti-neutrinos, cr' = (3cr /8)(2700£[5])/(77r 4 ) « (3/2)tr , and 
cr ~ 1.7 x 10~ 44 . Because both electron neutrinos are treated us- 
ing a single opacity, C A = 1/2 is set as applicable for the electron 
neutrino because their emission rate is generally larger or equal to 
the electron anti-neutrino emission rate when the neutrino chem- 
ical potential is zero. The electron-neutrino scattering expression 
above assumes the neutrino energy is thermal, which is inaccurate 
in the optically thin regime. However, this only leads to order unity 
corrections in our results for marginally optically thick regime. 
The absorption optical depth is 



(B92) 



(B85) 



where baryon reactions and photon-pair reactions have 
been collected together, and from Kirchhoff's law <r^ Ne = 

n+n+v+v— *n+n Q(i,n+n^>n+n+v+v 
°" v+v-te+H-er ~ Q»,e*+e-^v+vl '( n v cu 0,v), °~ v+v^y = 2o,r->v+v/( w vCH(),v)> 

where u Qv = (7/S)a r T 4 that includes both neutrinos and anti- 
neutrinos. Because the emission rate that is used merges both 
electron and pair capture processes, a single n b is used instead of 
each n p and n„ (for Y e = 1/2 this is accurate). Notice that n v ,n bf 
cancel out when obtaining the optical depth. 
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The total neutrino optical depth is T Vjt0 , = T v>sca + r y abs . As for 
photons, the two-stream approximation is used to obtain the proper- 
ties of the neutrinos for general optical depths. The p and r neutri- 
nos have been neglected up to this point because their optically thin 
emission rates are smaller than the electron types. However, in the 
optically thick regime all neutrino species have the same densities 
and pressure when assuming zero neutrino chemical potentials. To 
approximately capture this simple thermalization effect, note that 
the absorption and scattering opacities in the optically thick limit 
are similar to within an order of magnitude among neutrino species. 
So, the optical depth factors g[r v ] and h[r v ] are used to interpolate 
to the optically thick regime in order to include the p and t neutri- 
nos. The neutrino energy density loss rate is 

Gv = 7T-, ( B93 > 

where F v = cu , v h[T v ](l + 2/?[t„]), #„, abs = T y , abs /(dr y-llbs /d.s'), 

dT v ^bs/ds = ^bf^-tNe + ^"n+nmv-tn+n) ^v( ( - r V+v-*e + +e~ ^v+v^f) = 

Qo, v /(cu 0tV ). The (1 + 2A[t„]) factor in F v approximately accounts 
for p and t neutrinos in the optically thick regime by modifying 
the value of u v to be three times larger as required. The pressure, 
internal energy density, and number density of neutrinos are 

n v = « ,vg[r v ](l + 2g[T v ]\ (B94) 

U y = Ko,vg[T v ](l + 2g[T v ]), (B95) 

Pv = (1/3K, (B96) 

where n _ v = (45f [3]/(27r 4 ))(a y /k B )T 3 * (l/3)(a y /k B )T 3 . The 
(1 + 2g[r v ]) factor approximately accounts for p and t neutrinos 
in the optically thick regime. Future work can consider all neutrino 
species in detail and consider QED effects on the neutrinos that 
may be important when t„ > 1 . 



APPENDIX C: COLLISIONAL AND COLLISIONLESS 
RECONNECTION 

The goal of the following sections is to determine the thickness of 
the collisional current layer. The key point is that collisional recon- 
nection dominates unless dissipation (occurring on the scale of the 
layer thickness) is dominated by collisionless effects that occur on 
the scale of the plasma skin depth. 

The resistivity (77) is computed (section CI) in order to deter- 
mine the thickness of the collisional current layer (sections C2,C3). 
These calculations assume the background jet's values of b 2 and p h , 
the current sheet length L , and sheet separation length A . Then, 
collisionless reconnection is discussed (section C4). 

CI Classical Resistivity 

Consider the collisional resistivity on current-carrying electrons 
and positrons (with number density n e > x n em ) as due to: photon 
drag, proton collisions, or electron-positron collisions 8 9 . 

8 For the regimes of interest in this paper, the generalized Ohm's law shows 
that the current rise time for pairs (( r j se ~ tp/tg, the square of the elec- 
tron plasma time over the light crossing time of the current layer) is always 
shorter than the pair annihilation time (t a ~ (n em cr^ s cy 1 ) such that pairs al- 
ways contribute to the current. Protons are assumed to contribute negligibly 
to the current density, which is accurate for the regimes considered. 

9 The collisional resistivity computed here assumes all particles are non- 
degenerate, but the degeneracy of nucleons and electron-positrons is com- 



The resistivity is presumed to be associated with comov- 
ing 4-current j = en e iy d v d with comoving relative electron drift 
4-speed y d v d and electron charge e. Then, the resistivity r\' = 
(v ec m e )/(n e >e 2 ) can be determined by using the electron drift mo- 
mentum p e = y d v d m e , comoving electric field strength \E m \, ef- 
fective collisional frequency v ec , interaction force balance e|.E C0 | « 
PeV ec , and Ohm's law ej = e\E co \/rf. The magnetic diffusivity 
is then given by 77 = 77' c 2 /(4k), which can also be written as 
77 = d 2 ,v ec , where d e > = c/aj pe > is the skin depth of current- 
carrying electrons and positrons, oi p( ,< >non _ re i = ^4nn l ,ie 2 /m e is 
the current-carrying electron-positron non-relativistic plasma fre- 
quency, co pe , = 0V jn0n _ re iG^ (p 2 /4\ 2 _ 

1/2,1,3/2) /(2K 2 (/J)) is the asso- 
ciated relativistic plasma frequency, p = m e c 2 /(k%T e ), G() is the 
MeijerG function, and K2O is the modified Bessel function of the 
second kind (Bergman & Eliasson 2001). 

In addition to collisions between charged particles, current- 
carrying electrons and positrons with drift momentum p e = 
y d m e v d experience a photon radiative drag force of F e = 
-(4/3)(y d v d /c)u y <T S!i , where F e = e|£col- This gives an effective 
collisional frequency of v ec w F e /p e . From the magnetic diffusivity 
of 77 = d 2 ,v ec one obtains the resistivity due to photon drag of 

j) y ~ (4/3)d 2 (u y cr s]i c)/(m e c 2 ), (CI) 

(Goodman & Uzdensky 2008). 

The Coulomb resistivity for current-carrying electrons and 
positrons interacting with protons, electrons, or positrons is deter- 
mined by the mean free path given by A m f p ss l/((n pa iis + n p )cr), 
where in general rc pairs + n p > n p = n e so that pairs contribute 
an extra opacity that sets a lower-limit on the mean free path. The 
use of Mp a i rs + n p estimates the fact that electron-positron pair bath 
and the ion bath contribute (to order unity) the same to the per- 
particle Coulomb collision on the current-carrying electrons and 
positrons. For an effective collisional rate as given above (v ec ), the 
magnetic diffusivity is then 77 = ((rc pairs + n l ,)/n e r)(v d crm l ,c 2 )/(e 2 4T[) 
if each scatter is effective. If n pails — > or n e ' ~ 77 pairs » n p , 
then 77 = (v d crm e c 2 )/(e 2 4n), which is independent of the density 
of current-carrying electrons and positrons. In order to determine 
the resistivity, the cross section (cr) must be determined. 

Coulomb collisions have naive cross section of tr = nA 2 , 
where A c ss ^5/(l6Tc)e 2 /KE is the length scale over which 
Coulomb forces are important, where the kinetic energy of elec- 
trons in the center-of-momentum frame is given by KE ss 
^/( m e y th v th ) 2 c 2 + m 2 c A - m e c 2 , where v e ~ = v tb w 

c V©f'(2 + © c ')/(l + ® e ') is the thermal speed of electrons, y th = 
1/ Vl _ ( v th/c) 2 , and f - « (ujn e )/(m e c 2 ). There are additional 
relativistic corrections to A c of order unity not considered here (see, 
e.g., McKinley & Feshbach 1948). Also note that there are QED 
corrections that suppress the electron-electron cross section, but 
there is only weak suppression of electron-proton interactions that 
then dominate the Coulomb cross section (Storey & Melrose 1987; 
Sadooghi & Jalili 2007). For weakly coupled plasmas the Debye 
screening gives a corrected Coulomb cross section of 

a c x cr In A, (C2) 

puted and the fully degenerate electron-proton resistivity code by Potekhin 
et al. (1999) was used to check how electron-proton collisions are affected. 
Our work's results end up not depending upon electron-proton or pair- 
proton collisions, so the degeneracy effects (such as discussed in Rossi et al. 
2008 for electron-proton collisions) can be neglected in this study. Comput- 
ing radiative drag and pair drag on the degenerate current-carrying electrons 
and positrons is left for future work. 
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to an elongated reconnection layer. Plasma is forced into a confined 
space before being ejected out along a geometrically narrow chan- 
nel, which causes the reconnection rate to be slow 10 . 

The Petschek model uses the fact that topological changes in 
the electromagnetic field do not have to take place in the same re- 
gion where most of dissipation occurs. The Petschek model has a 
very small inner Sweet-Parker-like layer (called the central diffu- 
sion region) at the X-point where topological changes occur, but 
the current layer generates slow-mode shocks that dissipate most 
of the electromagnetic energy and deflect the incoming plasma into 
two relatively wide exhaust channels. The ejected plasma there- 
fore has a large geometrical exit channel allowing faster ingoing 
plasma and so faster reconnection. Petschek reconnection occurs 
with an asymptotic reconnection speed that typically relaxes to 
V r/ V A ~ 1/lnS (including in the relativistic regime; Lyubarsky 
2005). For a Lundquist number S » 1, the Petschek reconnec- 
tion rate is v r /vA ~ n/ (4ln(S (v r / va) 2 )), giving v r ~ 0.018c for 
5 = 10 22 for a highly magnetized jet (Lyubarsky 2005). 



Figure CI. Structure of Sweet-Parker and Petschek reconnection layers of 
length 2Lq: Top panel shows Sweet-Parker reconnection with narrow ejec- 
tion channel of thickness 2i5sp. Bottom panel shows Petschek reconnection 
with broad ejection channel between two shocks. The reconnection velocity 
is denoted by v r , and the ejection velocity by u (u ~ »a. the Alfven speed). 
In Petschek reconnection, most of the energy dissipation takes place in the 
slow-mode shocks, while in Sweet-Parker reconnection dissipation occurs 
within the central dissipation region. Both velocities are shown as green 
lines. The magnetic field, B and Bo, is shown as blue lines. Sweet-Parker- 
like reconnection occurs in collisional plasmas and is slow, while Petschek- 
like reconnection occurs in collisionless plasmas and is fast. 

where In A ss ln(A D /A c ) is the Coulomb logarithm and A D as 
Vth/ujpn' is the Debye length. 

The classical resistivity for current-carrying electrons and 
positrons interacting with protons, electrons, and positrons that 
accounts for Debye screening gives the Spitzer resistivity of 
r] s as dj,v c . While the relativistic expression is used, in the non- 
relativistic limit this becomes simply 

ns *[lS) CrM) ~ 3/2 ' nA ' (C3) 

where 6> f = (k^T) I (mgC 1 ) is the dimensionless electron temperature, 
and r e = e 2 /(m E c 2 ) is the classical electron radius. 
The total resistivity is then taken to be 

*1 = Vy + (C4) 

which is a simple sum because the resistivities act on the same 
current-carrying electrons and positrons. In most of parameter 
space for the GRB jets considered, photon drag dominates the 
Spitzer resistivity. Note that neutrino drag is generally negligible. 

C2 Reconnection Models 

Given the length (L ) of the current layer (as determined in 
section 5.1), one can derive two classes of solutions for mag- 
netic dissipation: the slow Sweet-Parker type reconnection (Sweet 
1958; Parker 1957, 1963a) and the fast Petschek-type reconnection 
(Petschek 1964). Which of these two modes of reconnection is ul- 
timately realized will be determined by the plasma collisionality as 
discussed in section 5.3. 

As shown in Figure CI, the Sweet-Parker regime corresponds 



C3 Collisional Reconnection 

Collisional reconnection is assumed to be approximately like 
the Sweet-Parker model of a stationary dissipating current sheet. 
Lyubarsky (2005) found that the relativistic Sweet-Parker recon- 
nection rate, layer thickness, and ejection velocity all behave as ex- 
pected from non-relativistic Sweet-Parker theory. For Sweet-Parker 
layer extent 2L , the Sweet-Parker layer has thickness 

6 SP ~ ,[^ = L S- il2 , (C5) 

which is associated with a reconnection velocity of 

v r , sp ~ y- =VaS- [/2 , (C6) 
OSP 

and an ejection velocity of u ~ v A . where recall that 5 = L v a /tj is 
the Lundquist number. The reconnection timescale is T sp ~ t a Vs, 
where ta = Lq/va is the fluid transit flow time along the length 
of the layer. For very large S , the SP solution may be unstable to 
plasmoid formation (see., e.g., Uzdensky et al. 2010). 

Uzdensky & McKinney (201 1) found that one can readily ex- 
tend the non-relativistic Sweet-Parker theory to include radiative 
cooling. The results in their section 2 apply for arbitrary optical 
depths. They obtained such results by using the energy equation 
rather than the incompressible assumption as done in Lyubarsky 
(2005). An interesting point is that the effects of cooling decouple 
from rest of equations, which allows one to obtain the Sweet-Parker 
type result from pressure equilibrium as usual and then indepen- 
dently obtain the constraint from radiative cooling conditions that 
determine a baryon compression ratio, A = n c /no, i.e., the ratio of 
plasma density in the center of the layer, n c , to the plasma density 
in the upstream region, no. 

The solution that Uzdensky & McKinney (201 1) obtained for 
general optical depths ultimately involves two equations: pressure 
equilibrium across the layer and radiative energy balance across 

10 Many effects are neglected, such as secondary instabilities (Dahlburg 
et al. 2003, 2005), 2D tearing instability (Loureiro et al. 2007; laroschek & 
Hoshino 2009) (but a temperature-dependent resistivity, non-uniform flow 
along the layer, and pressure anisotropy can stabilize the tearing mode; 
Biskamp 1986a; Priest & Forbes 2000), interchange instability (Uchida 
& Sakurai 1977), resistive kink mode (Schnack & Killeen 1979), an 
anisotropic Spitzer resistivity, realistic thermal conduction, and viscosity. 



© 201 1 RAS, MNRAS 000, 1-35 



32 J. C. McKinney, & D. A. Uzdensky 



the layer. While a global jet+layers structure is not sought in the 
present paper, a global solution would show that radiation pressure, 
like the magnetic pressure, would have its own structure across the 
jet. This justifies balancing the total gas+radiation pressure of the 
complex against the electromagnetic pressure of the jet. 

To obtain the radiative Sweet-Parker solution, first one im- 
poses the pressure balance condition 

Pem ~ P g (T), (C7) 

where p EM = b 2 /(Sn) is the electromagnetic pressure. This equa- 
tion corresponds to a loss of magnetic pressure within the layer that 
is recovered by a balanced thermal gas pressure. 

Second, the energy equation determines the effects of radiative 
transport (Uzdensky & McKinney 201 1). The result is that one sets 
the compression ratio of the rest-mass density within the current 
layer of size £ S p' (given below) to be 

A-^l, (C8) 

where recall that Q g (T) is the gas energy volume loss rate, and 

g SP = b 2 l(AnL lv A ), (C9) 

where v A is the usual upstream Alfven velocity (Uzdensky & McK- 
inney 201 1). The quantity gsp is the energy density dissipation rate 
if the layer corresponds to the standard thin Sweet-Parker solution 
with thickness S sp . Note that A does not depend upon the resistiv- 
ity. Equation (C8) applies only in the strong cooling regime, i.e., 
if A » 1. If A < 1, then one sets A = 1 corresponding to weak 
cooling limit giving back the original non-radiative incompressible 
Sweet-Parker solution. This A factor is applied to all baryon densi- 
ties (rib) and baryon-associated electron densities (n e ), except when 
used in opacity integrals because total baryonic mass across the 
jet is conserved regardless of localized compression in each thin 
current sheet. For example, there are no changes to the opacity in- 
tegrals over baryonic-associated mass when performed for equa- 
tions (B27, B31, B89, B92) 11 . 

Once the energy density loss rate (Q g ) is comparable to (or 
larger than) the energy dissipation rate (<2sp), the layer must com- 
press and undergo more rapid reconnection in order to balance the 
radiative losses. The reconnection velocity and current layer thick- 
ness are, respectively, given by 

*V,sp' ~ v r , sp A 112 , (CIO) 

6 sr > ~ S SP A- 112 , (Cll) 

(Uzdensky & McKinney 2011). This is only an approximation of 
a so-far undeveloped fully relativistic radiation reconnection the- 
ory. However, if A = 1 , then the radiation only contributes to the 
pressure and internal energy and the original non-relativistic Sweet- 
Parker solution should be accurate even in the relativistic regime 
(Lyubarsky 2005). In most of parameter space considered in the 
present study, A ~ 1 is found to hold. This implies that the usual 
Sweet-Parker solution is generally valid to order unity. 

The time-rate of change of reconnected magnetic flux is 

dP Idz 

— = \-E z \ = \v r B y /c , (C12) 
cdt 1 ' 1 

where B y ~ \b\ (the comoving electromagnetic field strength 

1 1 Because the free nucleon fraction depends upon density and not column 
density, for a pre-collapse baryon density pb one should integrate to get 
r oc AX R u C [Apb]S + X ntlc [pb](R j — 8). However, because generally Rj/6 »• A 
is found, this is a negligible effect and r oc X nuc [pt,JRj is set as usual. 



in the jet). In the non-relativistic incompressible limit, this 
gives a time-rate of change of magnetic flux of dP z /(dzcdt) ~ 
\v A S~ il2 \b\/c\ <k v A \b\/c, which implies a slow reconnection rate 
compared to Alfven timescales. The reconnection timescale is 
given by the time needed to reconnect a finite amount of flux of or- 
der \b\A dz present in the volume between multiple current sheets. 
This gives r r ~ Qb\A dz)/(dPJdt) ~ (A /« A ) VS, which is much 
longer than it takes an Alfven wave to cross A . 

In the relativistic regime the inflow speed may reach v r ~ v A ~ 
c and significant Lorentz contraction might occur when there is 
a significant loss of energy through radiation or if the Lundquist 
number 5 is smaller than the comoving magnetization b 2 / (&np b c 2 ) 
(Lyutikov & Uzdensky 2003). However, when including both the 
energy and momentum equations (Lyubarsky 2005), this is found 
not to occur in the Sweet-Parker regime. This validates the as- 
sumptions used in this work. Relativistic reconnection is an active 
area of research (Watanabe & Yokoyama 2006; Komissarov et al. 
2007; Hesse & Zenitani 2007; Zenitani & Hoshino 2007; Zenitani 
& Hesse 2008; Zenitani et al. 2009; Tenbarge et al. 2010), so as- 
pects of this work may require modification. 

C4 Collisionless Reconnection 

Sweet-Parker reconnection is slow because plasma has to flow 
through a narrow channel. Microphysics (e.g. anomalous resistiv- 
ity or the Hall effect) alone can not enhance the reconnection rate if 
the current layer is preserved in the Sweet-Parker configuration. For 
example, as the current layer thickness S — > di, Hall effects domi- 
nate the resistivity. The Sweet-Parker analysis yields a reconnection 
speed of v r ~ (di/L )v A , which still gives v r <K v A simply because 
of the assumed Sweet-Parker geometry. However, if a transition to 
collisionless reconnection occurs as discussed in section 5.3, then 
the Sweet-Parker geometry can be disrupted into a Petschek-like 
geometry leading to fast reconnection. 

One candidate for fast collisionless reconnection is the Hall 
effect. In the Generalized Ohm's law the Hall term is j x B/(n e >ec) 
operating on the ion skin depth where electrons and ions decouple. 
The Hall effect involves a two-fluid laminar reconnection configu- 
ration with v r < 0Av A (Mandt et al. 1994; Shay et al. 1998; Birn 
et al. 2001; Bhattacharjee et al. 2001; Cassak et al. 2005; Yamada 
et al. 2006; Daughton et al. 2006). 

Another candidate is an anomalous resistivity, which can be 
due to plasma microinstabilities (Parker 1963b; Coroniti & Eviatar 
1977; Tsuda & Ugai 1977; Parker 1979; Syrovatskii 1981; Biskamp 
1986b; Taylor 1986; Parker 1988b; Scholer 1989; Masuda et al. 
1994; Kulsrud 1998; Erkaev et al. 2001; Kulsrud 2001; Biskamp & 
Schwarz 2001; Malyshkin et al. 2005; Melrose 1986; Hasegawa 
1975; Begelman & Chiueh 1988). Anomalous resistivity can be 
triggered, e.g., when the drift velocity (v d ) exceeds some critical 
velocity (v c ), such that v d = j/(en e >) > v c ~ u t h- This process 
can drive plasma instabilities developing microturbulence where 
scattering of electrons by waves enhances resistivity. As the layer's 
thickness S decreases down to the critical thickness 6 C = cBI(Anj e ) 
where j e = en^y^v^, anomalous resistivity turns on (this scale is 
usually comparable to di). Anomalous resistivity allows for higher 
resistivity than Spitzer, but it also enables the Petschek-like geom- 
etry (see, e.g., Kulsrud 2001; Biskamp & Schwarz 2001; Yamada 
et al. 2006; Zweibel & Yamada 2009). 

Overall, collisionless Petschek-type reconnection operates on 
the ion skin depth 

di = c/cj pi , (CI 3) 
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where co pi is the ion plasma frequency. The ejection velocity is the 
same as in the Sweet-Parker case with velocity order v A . 
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